Quantum magnetism on the Cairo pentagonal lattice 
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We present an extensive analytical and numerical study of the antiferromagnetic Heisenberg model on the 
Cairo pentagonal lattice, the dual of the Shastry-Sutherland lattice with a close realization in the S = 5/2 com- 
pound Bi2Fe40<). We consider a model with two exchange couplings suggested by the symmetry of the lattice, 
and investigate the nature of the ground state as a function of their ratio x and the spin S. After establishing the 
classical phase diagram we switch on quantum mechanics in a gradual way that highlights the different role of 
quantum fluctuations on the two inequivalent sites of the lattice. The most important findings for S = 1/2 in- 
clude: (i) a surprising interplay between a collinear and a four-sublattice orthogonal phase due to an underlying 
order-by-disorder mechanism at small x (related to an emergent J1-J2 effective model with J2 S> Ji), and (ii) 
a non-magnetic and possibly spin-nematic phase with d-wave symmetry at intermediate x. 
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PACS numbers: 75.10.-b, 75.10.Hk, 75.10Jm, 75.10.Kt 



I. INTRODUCTION 

Geometric frustration is at the heart of strong correlations in 
many models of quantum magnetism. 1-3 Given that the most 
elementary building block of frustration is the triangle, lat- 
tice spin models with triangular units have been the minimal 
candidate models for realizing novel phases of matter in frus- 
trated magnetism and have thus been explored widely over 
the years. Perhaps the most celebrated example in two dimen- 
sions is the kagome lattice antiferromagnet (AFM), an array 
of corner-sharing triangles, where frustration leads to an ex- 
tensive number of classical ground states (GS). 1 As a result, 
quantum fluctuations play a non-trivial role and may even fa- 
vor a spin-liquid state, as it is currently believed for the ex- 
treme quantum-mechanical S=l/2 case. 1 

Another elementary unit with built-in frustration is the pen- 
tagon. As in the case of the triangle, an AFM exchange inter- 
action on a single pentagon favors a coplanar classical state 4 
which is frustrated, in the sense that not all sides of the pen- 
tagon are fully satisfied. This state is a spiral with pitch angle 
47r/5 which, compared to the 120° angle in the single triangle 
case, might suggest a lower degree of frustration. 

However, there is a generic aspect of pentagonal lattice 
models which underlies a deeper degree of frustration and 
complexity. This is related to the fact that there is no Bra- 
vais lattice of pentagons and so, unlike the 120°-state which 
can be easily "tiled" over the triangular lattice, it is not a priori 
evident that the single pentagon minimum is a good starting 
point for the description of the global low-energy physics. In 
fact, even in the finite-size case of the dodecahedron, which is 
a uniform tiling of pentagons on the geometry of a sphere, the 
GS is not related to the single pentagon minimum. 4 

A related aspect of 2D pentagonal lattices is that they of- 
ten consist of two or more inequivalent sites and bonds (see 
below). Together with the above, one then expects non-trivial 
classical and quantum-mechanical phases in pentagon-based 
lattice models, and this is a direction in frustrated magnetism 
which is largely unexplored. 

As mentioned above, the 2D plane has the generic prop- 



erty that cannot be tiled using regular pentagons. One alterna- 
tive, which is realized in the pentagonal Penrose lattice, is to 
combine pentagons with other motifs in order to fill the void 
spaces. 5 Another is to use irregular pentagons which leads to 
the 14 pentagonal tesselations known so far. 6,7 Among them, 
there are two tilings that have attracted interest in frustrated 
magnetism. The first is what is called "the pentagonal lattice" 
in some literature (see e.g. Fig. 2 in Ref. [8]), and the sec- 
ond is the "Cairo pentagonal lattice" 9-14 which is the dual of 
the Shastry-Sutherland lattice (see Fig. 9 of Ref. [12]), and is 
shown in Fig. 1. 

The Cairo lattice is the main subject of this study. Its main 
features can be seen in Fig. 1 . First, there are two inequivalent 
sites with coordination numbers 3 and 4, and likewise there 
are two inequivalent bonds connecting 3 -fold with 3 -fold sites 
(thin blue) and 3-fold with 4-fold sites (thick black). In partic- 




FIG. 1: (Color online) The Heisenberg model on the Cairo pentago- 
nal lattice considered in this article. Solid thick and thin bonds stand 
for the J33 and J43 exchange couplings. The (green) rectangle in the 
middle denotes the 6-site unit cell of the Cairo lattice. The four-fold 
coordinated sites form a tilted square lattice (denoted by thin dashed 
lines) which is further divided into two square sublattices denoted by 
large open (red) and filled (blue) circles. 
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ular, the 4-fold sites form a square lattice which is represented 
by dashed lines in Fig. 1. The Cairo lattice has a square Bra- 
vais lattice with a unit cell (green rectangle in the middle of 
Fig. 1) of six sites, four of which are 4-fold and two are 3-fold 
coordinated. 

Among the previous studies on the Cairo lattice, we point 
out the Ising model study by Urumov 9 and Rojas et al. 10 , and 
that by Ralko 11 who studied a Hubbard model with hard-core 
bosons (equivalent to an XXZ model under a staggered mag- 
netic field). Our own study focuses on the Heisenberg model 
described by the Hamiltonian 

H = Jij Si ■ Sj , (1) 

where the sum runs over the nearest-neighbors (ij) of the 
Cairo lattice. As shown in Fig. 1, we consider two different 
exchange interactions (which is the minimal number imposed 
by the symmetry of the Cairo lattice), J33 and J43, and focus 
on the regime where both are antiferromagnetic. By tuning 
the ratio x = J43/J33 and the spin 5, we shall be able to 
drive the system through a number of phases some of which 
have strong quantum-mechanical origin. As we are going to 
show below, the rich physics of this model is intimately con- 
nected to the presence of two inequivalent sites and bonds in 
the lattice. 

Besides the general motivation outlined above, there is an- 
other theoretical motivation for looking at this specific lattice 
model. This goes back to the work by Raman et al. 12 who 
proposed a connection between Quantum Dimer models 15 
(QDM) and the so-called Klein models using a well controlled 
decoration procedure on a number of frustrated lattices. The 
Cairo lattice has most ingredients for resonating valence bond 
physics, 15 since it is non-bipartite and has large resonance 
loops with even-length. 12 And indeed, one of the main pre- 
dictions of that study is that the QDM on the Cairo pentago- 
nal lattice has an extended spin liquid GS. This leads to the 
question whether such a phase would survive if one includes 
dynamics out of the singlet manifold, and so a study of the 
5 = 1/2 Heisenberg model is a natural extension. 

In parallel, it turns out that the Cairo pentagonal lattice is 
not only of purely theoretical interest, since the magnetism of 
Bi2Fe4C<9 offers a somewhat close realization of this model 
with 5 = 5/2. 13 ' 14 This compound was originally synthesized 
in the seventies by Shamir et al., 16 but it has attracted recent 
interest since it is a common by-product in the synthesis of the 
well known multiferroic compound BiFeOg. In fact, Singh et 
al., 11 have shown that Bi2Fe40g also shares some magneto- 
electric properties. 

A magnetic characterization in single crystals of Bi2Fe40g 
has been given by Ressouche et al. 13 Despite the large Curie- 
Weiss temperature 9 ~ 1670 K, this material orders magnet- 
ically at much lower temperatures Tn — 238 K, which is 
the standard signature of magnetic frustration. The most non- 
trivial finding is the nature of the low-T phase: It is a copla- 
nar configuration, whereby the 4-fold coordinated Fe 3+ spins 
form four orthogonal sublattices, while the 3-fold coordinated 
spins bind antiferromagnetically with each other and in the 



direction of the local exchange field exerted by their neigh- 
boring 4-fold coordinated sites. 

As explained by Ressouche et al., Bi2Fe40g is not a 
perfect realization of the Cairo lattice model. The first rea- 
son is that each pentagonal unit of the a6-plane comprises 
seven physical spins, since there are two ferromagnetically 
(FM) coupled Fe 3+ spins residing at each 4-fold site. Sec- 
ondly, the minimal microscopic model description of this 
compound comprises three in-plane and two out-of -plane ex- 
change pathways. 13 Despite this, the classical configuration of 
Bi 2 Fe40g seems to survive in a much larger parameter space, 
since the same state appears also in our more symmetric ver- 
sion of the model. 



II. MAIN RESULTS AND ORGANIZATION OF THE 
ARTICLE 

Our target is the GS phase diagram of the model in the x- 
5 plane. To accomplish this goal and establish a fairly good 
understanding of the various phases of the model we begin 
with the purely classical limit (Sec. Ill) and then we switch on 
quantum-mechanics in a gradual way by using different levels 
of approximations and complexity. In doing so we shall also 
be able to highlight the different role of quantum fluctuations 
at the two inequivalent sites of the model. 

Our theoretical predictions for the phase diagram in the x-S 
plane are presented in Fig. 2. The classical (large-5) phase- 
diagram is shown at the top line of Fig. 2. It consists of 
three magnetic phases: (i) the orthogonal phase (x < ^/2) 
found for Bi2Fe40g, (ii) a collinear 1/3-ferrimagnetic phase 
(x > 2), and (iii) an intermediate < x < 2) mixed phase 
which combines both (i) and (ii). We should note here that the 
1/3-ferrimagnetic phase appears also for Ising spins 10 and for 
hard-core bosons. 1 1 

Our first step to include quantum fluctuations in the prob- 
lem is by a standard linear spin wave theory (Sec. IV). Our re- 
sults for the renormalized spin lengths indicate three regions 
(see second line of Fig. 2) with strong quantum fluctuations: 
(i) at the 3-fold sites for small x, (ii) at the 4-fold sites as 
we approach the classical transition between the orthogonal 
and the mixed phase (x = V~2), and (iii) at the 4-fold sites 
throughout the mixed phase. These are the regions where new 
competing phases might appear as we approach the extreme 
5 = 1/2 quantum limit. By contrast the 1/3-ferrimagnetic 
phase seems to survive quantum fluctuations. 

One may immediately realize that the strong fluctuations at 
the 3-fold sites at small x result from the following generic 
feature. At small x, the 3-fold sites prefer to bind into 
quantum-mechanical singlets which are very different from 
the product up-down configuration of the classical ansatz. 
This strong tendency to form singlets on the J33 dimers can 
be captured by a modified variational ansatz (Sec. V) whereby 
the 3-fold dimers are treated fully quantum-mechanically and 
the 4-fold sites still classically. 

A numerical minimization of this ansatz for 5=1/2 gives 
the phase diagram shown in the third line of Fig. 2. The first 
important result is that the orthogonal state is now stabilized 
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FIG. 2: The Heisenberg model on the Cairo pentagonal lattice. Left panel: Evolution of the phase diagram as we go from the classical (top) 
to the quantum limit (bottom). Quantum fluctuations are included gradually, in a way that highlights the different role of the two inequivalent 
sites iof the lattice. The regions indicated by "SF" in the second and fourth lines are the ones with strong semiclassical fluctuations. Right 
panel: The main magnetic phases appearing in the phase diagram. The mixed phase (not shown) interpolates between the orthogonal (a) 
and the 1/3-ferrimagnetic phase (c) as described in the text. In (a) and (b) we also show our labeling scheme for the corresponding bosonic 
operators (for each of the 6 sites of the unit cell enclosed by the green solid line) that appear in LSW theory. Note that in the collinear phase 
(b), the spins in half of the J33 bonds (e.g. the ones labeled by c; and e*) do not feel any exchange field from the neighboring 4-fold sites. 



up to £ = 2, i.e. the mixed state does not survive in this ansatz 
for low enough spins S. In addition, the present ansatz pre- 
dicts a finite staggered polarization on the J33-dimers, apart 
from the strong tendency to form singlets. This can be ex- 
plained by the fact that the local exchange field exerted from 
the 4-fold sites in their orthogonal configuration is staggered 
and can thus admix a finite triplet |in) amplitude with the sin- 
glet wavefunction as soon as x is finite. 

We may go one step further and include stronger quantum 
fluctuations by performing an expansion around the orthog- 
onal variational GS found for x < 2. The natural way to 
do this is to perform a standard Holstein-Primakoff expansion 
for the 4-fold spins and a bond-wave expansion 1819 for the 3- 
fold dimers. The corresponding quadratic theory is presented 
in Sec. VI and its predictions for S = 1/2 are shown in the 
fourth line of Fig. 2. The main result is the presence of very 
strong fluctuations above x ~ 1, which provides strong evi- 
dence that the variational treatment is not a good starting point 
in the window 1 < x < 2, and that a new phase will be stabi- 
lized in this regime for low enough spins S. 

The last place to look for quantum fluctuations is at the 4- 
fold sites (Sec. VII), which have been treated classically or 
semi-classically up to now. It turns out that these fluctua- 
tions are responsible for some very rich physics at low en- 
ergies. The main reason for this is that the x = limit has a 
highly degenerate GS manifold since the 4-fold sites are free 
to point up or down in this limit. As a result, the low-energy 
physics at small x is governed by effective interactions be- 
tween the 4-fold sites which are mediated by the virtual fluc- 
tuations of the J33 -dimers out of their singlet GS. By inte- 
grating out these high-energy fluctuations up to fourth order 
in x, we have derived an effective low-energy theory for spins 



S = 1/2, which governs the interactions between the 4-fold 
sites up to x ~ 0.4 — 0.6. 

The first key result of the effective theory is that the nearest- 
neighbor coupling 3\ is much smaller than the next-nearest- 
neighbor coupling J2 due to destructive quantum interfer- 
ences. Thus, at small enough x, the effective model reduces 
to the well studied 20-25 J\-Ji AFM model on the square lat- 
tice with J 2 Ji . An immediate consequence is that there 
must exist a critical value of x (see fifth line of Fig. 2) below 
which the orthogonal phase becomes unstable in favor of the 
collinear phase through an order-by-disorder mechanism. 

The second important result from the effective theory is the 
appearance, in fourth-order in x, of a 4-spin coupling term K 
which involves the 4 spins in every plaquette of the square lat- 
tice, and is similar to the well-known ring exchange term. 26 
Our exact diagonalizations provide evidence that this plaque- 
tte term is the one that actually drives the system into the inter- 
mediate (1 < x < 2) quantum-mechanical phase mentioned 
above. 

It is worth noting that the collinear phase cannot not be 
detected using the above linear semiclassical theories. The 
reason is that the collinear phase is not the GS of either the 
fully classical or the variational ansatz, and therefore one must 
include interactions between the spin-waves or spin+bond 
waves respectively in order to stabilize this phase. To demon- 
strate this, we have employed (Sec. VIII) a spin-wave ex- 
pansion around the collinear phase using the effective model 
Hamiltonian and keeping up to quartic terms in the interac- 
tions. The results from a self-consistent mean-field decou- 
pling show that the collinear phase can indeed be stabilized at 
low enough x (fifth line of Fig. 2). 

It is by now quite clear that the phase diagram becomes 
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richer as we approach the extreme S — 1/2 limit. We have 
seen for instance that out of the three classical phases, the 
mixed phase does not survive quantum fluctuations for low 
enough spins S. We have also uncovered an interesting order- 
by-disorder mechanism which is at play in the small-x regime 
and gives rise to the collinear phase for S = 1/2. We have 
also learned from the spin+bond-wave expansion that a new 
phase is to be expected at intermediate x for low enough S. 
However the nature of this phase is not known yet. In par- 
ticular, one important question is whether the small x limit is 
a good perturbative limit for the description of this phase or 
not. If the answer is yes we must next identify the effective 
term that drives the transition, which is a difficult task given 
that our fourth-order effective model should be valid only up 
to x ~ 0.4 - 0.6. 

To solve these open issues we look at the full quantum- 
mechanical S = 1/2 problem using exact diagonalizations 
(ED) on finite-size clusters (Sec. IX). Our numerical results 
on the original Cairo lattice model confirm the presence of 
two nearly decoupled AFM sublattices for small x which 
is the physics we expect from the effective model at large 
J2/J1. Both the symmetries of the low-energy spectra and 
the GS correlations show signatures of the orthogonal and the 
collinear phase, but we are not able to pinpoint the exact tran- 
sition between the two phases given that the locking between 
the two AFM sublattices present at large J2 / J\ takes place at 
very large length scales. 27 ' 28 

The results also establish the presence of a new phase at in- 
termediate values of x before we reach the 1/3-ferrimagnetic 
phase. Apart from a GS level crossing we also find a new 
spectral structure at low energies. In particular, the magneti- 
zation process in a field now shows steps of AS Z = 2 which is 
typical for collinear spin-nematic phases. 29,30 Moreover, this 
phase seems to be adiabatically connected to the GS manifold 
of the x = limit, which suggests that it is driven by one of 
the couplings in the effective model. 

To clarify this issue we have also performed ED in the 
effective model on the square lattice but with unconstrained 
J1/J2 an d K/J 2 - We have found that the low-energy spectral 
structure in the regime K 3> J12 is very similar to the one in 
the Cairo lattice model at intermediate x. This suggests that 
the intermediate phase is driven by one of the three topologi- 
cally different 4-spin exchange terms of the model, which we 
identify as the one which has an enlarged SU(2)xSU(2) sym- 
metry. We show that this term favors a spin-nematic phase 
with d-wave symmetry similar to the one found by Shannon 
et al. 29 The main difference is that here the symmetry break- 
ing seems to take place in one of the two sublattices of the 
square lattice only, while the symmetry of the other sublattice 
remains intact. 

The remaining part of the article is organized along the 
main lines described above. We shall mainly focus on the 
central aspects and predictions of each separate approach and 
relegate technical details to the Appendices. 



0.05 




. 1 1 I 1 . 1 . . 1 . I 1 1 1 . I . 1 . . 1 a I 

1 1.2 1.4 1,6 1.8 2 2.2 

X=J 43 /J 33 



FIG. 3: Classical GS energy per site (measured from the energy of 
the orthogonal phase) as a function of x — J43 / J33. The lines stand 
for the analytical expressions given in Eqs. (7-8), and the symbols 
are numerical data from CMC simulations at /3J33 = 10 3 . 



III. CLASSICAL LIMIT 

In the classical limit we find three different GS's as a func- 
tion of x. At large x the lattice becomes effectively bipartite 
and one can minimize the energy by a collinear arrangement 
of up and down spins on the 3-fold and the 4-fold sites respec- 
tively (see Fig. 2(c)). Since the number of 3-fold sites is twice 
the number of 4-fold sites, N 3 = 2/Y4, this is a ferrimagnetic 
configuration with a total magnetization of 1/3. This phase 
remains stable down to x = 2. 

In the opposite limit of small x, the classical GS is the or- 
thogonal configuration shown in Fig. 2(a), which is the one 
found for Bi 2 Fe40g. 13 Here the 4-fold sites form an orthog- 
onal 4-sublattice configuration, while the nearest neighbor 
3-fold sites align antiferromagnetically to each other and at 
the same time point opposite to the total exchange field ex- 
erted by the neighboring 4-fold sites. This phase remains sta- 
ble up to x = y/2. We note here that the orthogonal con- 
figuration of the 4-fold sites has been found previously on 
some ring-exchange 26 models on the square lattice, see e.g. 
Refs. [31,32]. 

This leaves a window between x = \/2 and 2 where the 
spins find a compromise between the two phases by combin- 
ing both into a mixed phase. Namely, at x = \/2 the 3-fold 
and 4-fold sites begin to tilt out of the plane but in opposite di- 
rections to each other. In particular, the projection of this non- 
coplanar configuration onto the xy-plane gives the orthogonal 
phase while the projection along the z-axis gives the 1/3 fer- 
rimagnetic state. At x = 2 the spins are completely aligned 
along the z-axis. 

All three classical phases are special cases of the ansatz (see 
right panel of Fig. 2(a)) 

S Qi = S {pi cos 9 e y — sin 9 e z ) (2) 

S/ 4 = S (piCOs9 e x — sin 6* e 2 ) (3) 

S bz = S di = S (-Pi cos 9' e+ + sin0' e 2 ) (4) 

S C! =S Ei = S {pi cos 9' e_ +sm6'e z ) (5) 
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The two 



where ^ = e l( ^' Ri , Q = (tt,tt), and e± 
angles 6 and 6*' account for the tilting out of the xy -plane of the 
4-fold and the 3 -fold sites respectively. For the orthogonal and 
the 1/3-ferrimagnetic phases, 6 = 8' = and 6 = 6' = tt/2 
respectively, while for the mixed phase sin 6 = y/2 — 4/x 2 , 
and sin 6' = § sin 6. The corresponding energies per site are 



Eortho/S 2 = -(1+2V220/3 
Smixed/ S = — 1 — X /3 , 
Eferri/S 2 = (1 - 4z)/3 . 



(6) 
(7) 
(8) 



To confirm that these phases correspond to the global min- 
ima we have also performed Classical Monte Carlo (CMC) 
calculations at low temperatures using the Metropolis algo- 
rithm. The average energies per site in some representative 
low-T equilibrium ensembles are shown in Fig. 3 for a num- 
ber of cluster sizes. The results are in excellent agreement 
with the above picture. In particular they confirm the pres- 
ence of the intermediate mixed phase. 

We note here that the energy of the mixed phase is only 
slightly below the energies of the neighboring orthogonal and 
ferrimagnetic phases, and this suggests that this phase may be 
quite fragile against quantum fluctuations. Indeed, the varia- 
tional treatment presented below in Sec. V will demonstrate 
that the mixed phase is unstable for S — 1/2 as soon as we 
include the quantum fluctuations on the 3 -fold sites. 



IV. SEMICLASSICAL EXPANSION 

A. Linear Spin-Wave theory 

Our next step is to assess the strength of quantum fluctua- 
tions by performing a separate semiclassical expansion around 
each of the three classical GS's. The Cairo pentagonal lat- 
tice has 6 sites (two 4-fold and four 3 -fold coordinated) per 
unit cell. So we introduce six bosonic operators denoted by 
ai,bi,Ci,di,ei, fi to describe the transverse fluctuations on 
each site of the i-th unit cell. For each site we define a local 
quantization z-axis (see Fig. l(a)-(c)) and perform a standard 
Holstein-Primakoff expansion, namely 



S 



2S cu 



(9) 



and similarly for the remaining sites of the unit cell. The 
standard procedure 33,34 for the diagonalization of the resulting 
quadratic Hamiltonian is outlined in App. A. Here we shall fo- 
cus on the two main quantities of interest: the renormalization 
of the GS energy and the renormalization of the local spin 
lengths by harmonic fluctuations. 



B. Results 

In the following we present the predictions of LSW theory 
for spins S = 1/2. We first discuss the harmonic correc- 
tion SE to the GS energy. This is shown in the upper panel of 



Fig. 4, where we also make a comparison to exact diagonaliza- 
tion results on clusters with 12 and 24 sites and S=l/2 (open 
symbols). We see that by including harmonic corrections to 
the classical energy one obtains a remarkably good agreement 
with the full quantum-mechanical results on finite-size clus- 
ters (the agreement becomes better with increasing number of 
sites especially at larger x). However as we are going to show 
below this agreement is quite deceptive for S = 1/2. Al- 
though the semiclassical theory accounts for most of the en- 
ergy, there is quite a lot of rich physics taking place at much 
smaller energies compared to the bare scale J 33 , and this hap- 
pens because the x = limit is singular, in the sense that it 
has a highly-degenerate GS manifold (see Sec. VII below). 

We next examine the quadratic correction SS to the local 
spin lengths of the two inequivalent sites of the lattice, which 
are shown in the lower panel of Fig. 4. The first important 
feature to note is the upturn of SS for the 3-fold sites as we ap- 
proach the limit x = 0. This is a manifestation of the singular 
nature of the x = limit mentioned above. In the semiclas- 



ferrimagnetic 




FIG. 4: Upper panel: Classical (E c i) and semiclassical (E c i + SE) 
GS energy per site from LSW theory around the orthogonal (black), 
the mixed (blue), and the ferrimagnetic state (red), as a function of 
x = J43/J33. For comparison, we also show (with symbols) the 
GS energy per site obtained from a fully quantum-mechanical nu- 
merical calculation (ED) on two clusters with 12 and 24 sites with 
periodic boundary conditions. Lower panel: Correction to the local 
spin length of the 4-fold (circles) and the 3-fold (squares) coordi- 
nated sites. The (green) diamonds show the results from LSW theory 
around the orthogonal phase in the effective model (see Sec. VII). 
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sical treatment we are expanding around the product up-down 
state on the J 33 -dimers, and thus we cannot capture the actual 
tendency to form quantum-mechanical singlet wavefunctions 
at small x. Below we shall correct for this effect for the case 
of S = 1/2 by integrating out the large J33 energy scale of the 
,/33-singlets and by deriving an effective Hamiltonian model 
for the 4-fold sites only. In doing this it will become imme- 
diately apparent that the orthogonal phase is actually unstable 
at small x and S — 1/2 against a collinear magnetic phase 
which is stabilized by an underlying order-by-disorder effect. 

A second feature in our results for 5S is its upturn for both 
types of sites around the transition between the orthogonal and 
the mixed phase. This is clearly a sign of strong quantum 
fluctuations, and suggests that another possibly non- magnetic 
phase might be stabilized for S = 1/2 in this regime. Our ED 
results for S — 1/2 will indeed provide strong evidence for a 
non-magnetic state in this regime. In particular, as we show in 
the following section, the mixed phase becomes unstable al- 
together by including quantum fluctuations on the J33 dimers 
only. 

Third, the harmonic corrections in the ferrimagnetic phase 
seem to be almost independent of x and suggest that this phase 
probably survives quantum fluctuations. Indeed our ED re- 
sults confirm this. 

Finally we should note that, except for the special region 
around x = 0, the correction for the 4-fold sites is generally 
larger than that of the 3-fold sites. For x > V2 they differ 
by about a factor of two. This is a manifestation of the very 
different role of quantum fluctuations on the two inequivalent 
sites of the lattice. It is also at odds with the simple intuition 
that higher coordination sites tend to behave more classically, 
but seems to be a consistent feature in a number of 2D lattices 
(see e.g. Ref. [35]). 



V. SWITCHING ON QUANTUM MECHANICS ON THE 
3-FOLD DIMERS: VARIATIONAL ANSATZ 

In the purely classical description one treats each spin as 
a classical vector pointing in some fixed direction in spin 
space. However, as we discussed above, in the limit of small 
x, the nearest neighbor 3-fold sites prefer to bind into sin- 
glets, which are locally entangled quantum-mechanical states 
that are very different from the up-down (or down-up) product 
state of the classical ansatz. In particular, the latter contributes 
an energy of — S 2 , while the singlet wavefunction contributes 
—S(S + l). 61 In order to capture this strong-singlet physics at 
small x we introduce a variational wavefunction |^ V ar) which 
still treats the 4-fold sites as classical vectors but leaves com- 
plete quantum-mechanical freedom for the 3-fold sites. Fur- 
thermore, to include both the orthogonal and the mixed phase 
as special cases of this ansatz we assume that the unit cell of 
l^var) is twice the unit cell of the Hamiltonian with ordering 
wavevector at (tt, it). In addition, we do not put any restric- 
tion on the directions of the 4-fold classical vectors. Thus 
the variational parameters in this ansatz are the 6 polar an- 
gles corresponding to three 4-fold classical vectors (by global 
SO(3) symmetry, the fourth 4-fold vector is forced to point 
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FIG. 5: Comparison between the numerical minimization of the vari- 
ational ansatz described in the text (symbols) and the analytical pre- 
diction (solid lines) from the problem of an AFM dimer in the pres- 
ence of a staggered field h s = \/2Sx (see App. B). Circles show the 
squared overlap between the optimal quantum-mechanical state | $) 
of a dimer and the full singlet state, and squares show the local spin 
lengths on the J33-bonds. 



in a fixed direction). In effect, this ansatz amounts to solving 
quantum-mechanically, for each set of the above parameters, 
the problem of an AFM S=l/2 dimer in the presence of two 
local fields of arbitrary directions. 

The first major result from the variational minimization is 
that the optimal GS corresponds to the orthogonal coplanar 
configuration up to 1 = 2 and to the ferrimagnetic state 
(where the J33 bonds form |ft) triplets) for x > 2. Hence 
including quantum-mechanics on the 3-fold sites makes the 
mixed phase unstable for low enough spin S, which confirms 
our expectation that this compromise between the orthogonal 
and the ferrimagnetic phase is fragile. 

We have double-checked this important result by search- 
ing for the minimum energy of the same variational problem 
but in the more restricted parameter space whereby the 4-fold 
classical spins are tilted away from the ly-plane by an an- 
gle 8 (i.e. as in the mixed phase configuration). This angle 
is the only variational parameter which makes the problem 
of finding the minimum much more tractable numerically. In 
fact, this problem is equivalent to that of an AFM dimer in 
the presence of a staggered field h s = \J~2xS cos 8 e x , plus a 
uniform field h„ = 2xSsin8 e z , i.e. it is an extension of the 
staggered field-only case treated in App. B. And indeed, the 
GS of this problem has 8 — (coplanar phase) for < x < 2, 
but 8 = 7t/2 (ferrimagnetic phase) for x > 2. 

We now look at the main quantities of interest in the above 
optimal variational state for S = 1/2. The first is the overlap 
of the optimal GS on the dimer with the exact singlet wave- 
function. The second is the polarization of the two 3-fold sites 
forming a dimer. Both quantities signify the amount of triplet 
admixture in the GS, and they are shown in Fig. 5, where they 
are also compared to the corresponding analytical predictions 
given in App. B. Our first comment is the overlap remains 
quite large up to x = 2 which signifies that the strong x = 
coupling limit is a good perturbative limit for discussing the 
physics of the full S — 1/2 quantum problem even at interme- 
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FIG. 6: Bosonic operators and local quantization axes (for the four- 
fold sites and the first site of each dimer) in the orthogonal state, 
where all spins lie on the xz-plane. 



diate x. Another feature is that the polarization of the 3-fold 
sites becomes immediately finite as soon as we switch on a fi- 
nite x. This happens because we are dealing with a quantum- 
mechanical dimer in the presence of a staggered field (i.e., the 
local exchange fields at the two sites of the dimer are oppo- 
site to each other in the coplanar phase) which admixes the 
triplet |£rj) into the singlet GS as soon as x is finite (see details 
in App. B). This is in contrast to the case of a AFM dimer 
in a uniform field which can polarize the dimer only above a 
critical value which is set by the singlet-triplet gap. 



VI. QUADRATIC FLUCTUATIONS AROUND THE 
VARIATIONAL ANSATZ 

A. Linear Spin+Bond-Wave Theory 

In the previous variational wavefunction the J33-dimers 
were treated fully quantum-mechanically but the 4-fold sites 
were treated as classical vectors. So we refine our treatment to 
include the quadratic fluctuations around the variational state, 
by performing a semi-classical spin wave expansion for the 4- 
fold sites and a bond-wave expansion 18,19 for the J33-dimers. 

Figure 6 shows the unit cell of the lattice and the orthogonal 
variational state around which we expand. For the 4-fold sites 
a,i and fi, we perform a Holstein-Primakoff expansion S*. = 



S — a~[a,i, — v25 ca (and similarly for S/ ; ) where we 
use the local quantization axes shown in Fig. 6. 

Using the above representation for the 4-fold spins, we may 
regroup the various terms in the Hamiltonian as follows 



H =H {33) +n\ 



(43) 



H 



(43) 



Ho + H ( 2 43) (10) 



where "H*- 33 - 1 contains the J33 coupling terms, while "H^ 43 -* 
contains the parts that come from the constant S from the z- 
componentof the 4-fold spins, and contains the remain- 
ing portion from the J43 coupling terms. More explicitly, 



Ho 



(S^i — S2,i) 



^33 S' M • S' 2ti - h s z s ^i ■ [S[ ti - S' 2 i 



(11) 



where h s = y/2xS. Therefore T-L describes two independent 
dimers in the presence of staggered fields, a problem that is 
solved in App. B. 

To account for the fluctuations driven by the remaining por- 
tion 'H 2 43 ' ) °f tne Hamiltonian we shall make a bond-wave 
expansion. We first introduce eight bosonic operators which 
create the singlet |s) = 4= (|t|) - | It)) and the three triplet 

states |*0 = |tt), \t-i) = III) and \t ) = -± (|||) + |||)), 
for each of the two J33 -dimers per unit cell. These bosons will 
be denoted by {s, t\,to, t^i}i and {s' , t^,^, (see Fig. 
6). The spin operators of the two sites of each dimer have the 
following bosonic representation 18 ' 19 

5f )2 = ii^s + s+toJ+i^J-ti-tlit-i) (12) 

5+ 2 = (4s - s+t-r) + ±= (t+to + i + t-i)C 13 ) 

We also introduce the bosons which create the two lowest 
eigenstates of Ho (see Eqs. (B2) and (B3)) for the two sep- 
arate dimers in the unit cell, namely 

V> + = iis + + d( c + , ip£ = v s + - u i + , (14) 
ip[ + = us'++vt' Q + , = v s'+ - u t'+ , (15) 

where the constants u,v are defined in App. B. To expand 
around the variational GS we take the constraint V'fV'i + 
tpt^z + + t-it-i = 1> and replace the right hand side 
with a large number M. We then assume that the bosons ipi 
and ip[ are condensed and perform a large-M expansion 19 



1p! ~Vm- 



1 



2VM 



W> + V>2+* + ii+*ii*-i) 



(16) 



Replacing in the Hamiltonian and keeping only quadratic 
terms in the expansion we arrive at 



n-E 



9 A k 



D+ 



(17) 



where E /N uc = (2M + 3)ei - e 2 - 1/2 - AxvuM/y/2, 
N uc = 6N is the number of unit cells, 61,2 are the single- 
dimer energies given in Eqs. (B2) and (B3), 
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FIG. 7: Main results from the spin+bond-wave calculation which captures the quadratic fluctuations around the variational GS. Left panel: The 
8 branches of hybridized spin+bond wave excitations at x — 0.5 (see text). Middle panel: GS energy per site as a function of x. Comparison 
between linear spin wave theories around the orthogonal and the mixed state (dashed lines), the spin+bond-wave theory, and ED data from 
the 12 and 24 site clusters. Right panel: Various GS expectation values as a function of x. stands for the local spin length of the 4-fold 
sites. (5i)o and (Si ) stand for the staggered polarization on the ,/33-dimers in the variational GS and when we include quadratic fluctuations, 
respectively, and similarly for (Si -82)0 and (Si 'S2). The symbols are ED results from the 12 (triangles) and 24 (circles) site clusters. 



= ( at, fi, Vai*, tl, k : *ti,fc> ^atfc' a - fc ' /- fe ' ^2,-*, ti,-*, t-i,-*, ^ 2 ,-fc> f-i,-*, ) > 

, (18) 

I 



and the 8x8 matrices Cfc and are given explicitly in 
App. C. To diagonalize this quadratic Hamiltonian we search 
for a Bogoliubov transformation = Vfe • A& as described 
in App. A, in terms of new bosons //-, ^2> • • • , u for 
which 

n ~ Eq + - J2 + • ■ • + ^sk) 
k 

+ J2 ("ik ^ a k + ... + uj 8 k ■ (19) 

k 

B. Results 

In the left panel of Fig. 7 we show the eight spin+bond- 
wave branches of excitations along some symmetry directions 
in the BZ and for x = 0.5. The spectrum consists of 2 
low-lying modes and 6 high-energy modes around E = J33. 
These modes arise respectively from the two spin-wave modes 
of the 4-fold sites and the 6 J33 -triplet modes per each unit 
cell. At x = 0, both the 4-fold sites and the J 33 -dimers are 
isolated and thus all modes are completely localized. In par- 
ticular, the two spin-wave modes have zero energy and the 6 
triplet excitations cost energy J33. The hybridization caused 
by a finite x then gives rise to the dispersion structure shown 
in the left panel of Fig. 7. 

Of particular interest are the Goldstone modes which ap- 
pear in the spectrum due to the fact that the variational state 
breaks the continuous SO(3) symmetry. From the nature of 
the variational state one expects two Goldstone modes, one at 
k = and another at the ordering wave-vector Q = (71", tt). 
However, the left panel of Fig. 7 shows that we actually have 



3 gapless modes, one at zero momentum and two at Q. The 
extra gapless mode at Q appears also in the linear spin-wave 
dispersions (not shown here) around the orthogonal state and 
is spurious in both cases. Such a spurious gapless mode ap- 
pears also in the J1-J2 model around the collinear phase in the 
large J2 regime where it may be lifted by including nonlinear 
terms. 36 In the present case, these non-linear terms actually 
select a different state than the orthogonal one. 
We now turn to the quadratic GS energy given by 

E = E ° + 2 ^ + ■■■+ W8fc ) ' (20) 

A; 

Our results as a function of x are shown in the middle panel 
of Fig. 7, where we also make a comparison with the semi- 
classical spin-wave energies from Sec. IV. Comparing with 
ED data on finite clusters with 12 and 24 sites, we see that the 
present Spin+Bond-Wave theory gives a better agreement than 
the Spin- Wave theory especially at small x. However at larger 
values of x, the mixed spin+bond-wave expansion does not 
deliver a better energy than the pure semiclassical expansion. 

Let us now look at the spin length of the 4-fold sites which 
is given by 

1 16 

uc k n=9 

and can also be calculated numerically by integrating over 
the BZ. The result is shown in the right panel of Fig. 7 as a 
function of x. We see that quadratic fluctuations destroy com- 
pletely the long-range order at the 4-fold sites around x ~ 1.5. 
This means that there exist quite severe quantum fluctuations 
in this regime which might destabilize the anstaz phase. 
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A similar conclusion arises by looking at the behavior of 
the exchange energy on the J 33 -dimers as well as their stag- 
gered polarization. These quantities can be calculated using 
the expressions 

(Si-S 2 ) = ^- u 2 + {2u 2 - l){i>+^ 2 ) +2u 2 {t+t 1 \22) 
(Sl 2 ) = ±uv (1 - 2(^2) - nth)) (23) 

where we have used (t^ti) = (t^t-i) (due to time reversal 
symmetry), and 

16 

(v> 2 + v> 2 ) = t-EEi^( 3 '")i 2 c 24 ) 

" uc k n=9 
1 16 

uc k n=9 

The results are shown in the right panel of Fig. 7 where they 
are compared to the corresponding behavior in the variational 
GS, i.e., without quantum fluctuations. We see that the ex- 
change energy becomes quickly reduced in magnitude with x, 
and it even crosses over to positive values above a; ~ 1.8. This 
behavior is drastically different from the corresponding result 
(Si • 82)0 in the variational GS and from the ED data. Thus 
the effect of quadratic fluctuations is quite strong. 

VII. FULL QUANTUM 5 = 1/2 LIMIT: EFFECTIVE 
LOW-ENERGY THEORY FOR SMALL x = J43 / J33 

Up to now, the spins at the 4-fold sites were treated classi- 
cally or semi-classically. However it turns out that there are 
quite strong quantum-mechanical effects on these sites which 
modify to a large extent the picture we have so far and brings 
rich physics which takes place at a smaller energy scale. 

At J43 = the system consists of isolated J33-dimers 
which form singlets, and free spins (4-fold sites) which are 
free to point up or down. This highly degenerate GS manifold 
is lower in energy from excited states by at least J33 which is 
the cost of promoting one singlet into a triplet. By switching 
on a small x the 4-fold spins begin to interact with each other 
through the virtual fluctuations of the 3 -fold dimers out of the 
singlet manifold. By integrating out these fluctuations one 
may derive an effective model description for the low energy 
sector using degenerate perturbation theory. Specific details 
of this expansion are provided in App. D. Here we summarize 
the main results. 

It turns out that an expansion up to fourth order in x pro- 
vides the essential low-E physics at small x. The major in- 
teractions between the 4-fold sites are depicted in Fig. 8. The 
first important insight comes already in second order pertur- 
bation theory. Namely that the effective n.n. coupling Ji van- 
ishes due to quantum interference since the two possible paths 
shown in Fig. 8 have opposite amplitude. The overall minus 
sign originates in the fact that the singlet wavefunction is an- 
tisymmetric with respect to the interchange of the two sites 
of the J33-dimer. By contrast there is only one possible path 




FIG. 8: Left: Lowest order processes in the effective model. The 
4-fold sites form a tilted square lattice (dashed lines). In second or- 
der in x, Ji vanishes because of quantum interference between two 
different paths (shown by arrows) with opposite amplitude. By con- 
trast, J2 is finite since there is a single path. In fourth order in x, a 4- 
spin term appears that invokes the 4 spins around a square plaquette. 
Right: Dependence of Ji, J2 and K on x, as given by Eqs. (28)-(30). 



connecting n.n.n sites and this gives a finite J 2 coupling in 
2nd order. A nonzero Ji appears first in third order together 
with a renormalization for J 2 . In fourth order, in addition to a 
renormalization of J\ and J 2 , one also obtains a four-spin ex- 
change term that involves the spins of each square plaquette. 
For a plaquette with a horizontal dimer (see the 7th-cluster in 
Table II), this term reads 

£ = (S r S 2 )(S3-S 4 )— (S r S4)(S 2 -S3)+(SrS 3 )(S 2 -S4) (26) 
or, in pictorial form, 




4 3434 3 



= ICh - -jfcv + &x • (27) 

This term is different from the usual ring exchange process 26 
on the square lattice. In particular, the coefficients of ICh and 
1C V are different here (the actual relative factor of —1/2 is not 
generic but is expected to be modified in higher orders of per- 
turbation theory), and this reflects the lack of the C\ symmetry 
around the center of each square plaquette in the underlying 
parent Hamiltonian on the Cairo lattice. In essence the full 
symmetry of the underlying Cairo lattice is only manifest in 
fourth order of perturbation theory, which is a good reason as 
to why one should push the perturbation theory at least up to 
4th order. Finally the coefficient of the last term IC X is the 
same as that of ICh and this reflects the underlying symmetry 
of exchanging S 2 and S3 (or Si and S4) in all 4-th order pro- 
cesses. However this is not a generic feature although it might 
still hold up to some order higher than fourth. 

Altogether, up to fourth order in x we have a low-energy 
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X = J 43 /J 33 



FIG. 9: Left panel: Self-consistent solution for the mean-field pairing fields that appear in the interacting spin wave theory around the collinear 
phase. Middle panel: The energy (in units of J2(x')) of the second excitation mode at k = and at k = (n, it). Inset: The BZ of the effective 
model is enclosed by the dashed (blue) lines and is set by the reciprocal vectors G x and G y . At small enough x the effective model reduces to 
the J1-J2 model on the square lattice whose unit cell contains one site only. Its BZ is enclosed by the solid (red) lines and is set by the vectors 
Gi and G2. Right panel: Spin-wave dispersions from NLSWT around the collinear phase for two representative values of x. 



effective Ji-J 2 -K model on the square lattice with parameters 

Ji = \ x3 + \ x4 (28) 
■h = \x 2 + ^ 3 - ^ 4 (29) 
K = x 4 (30) 

up to fourth order in x. As can be seen in the right panel of 
Fig. 8, these functions show a rapid increase above x ~ 0.4, 
so we expect that the fourth order series expansion does not 
converge beyond x ~ 0.4 — 0.6. 

From the above effective theory we get the following in- 
sights into the low-energy physics of the problem at small 
x. The first insight comes from the fact that we are dealing 
with a dominant J 2 coupling in the square lattice. This gives 
two square sublattices which are decoupled from each other 
and order antiferromagnetically. Now the angle 9 between the 
two Neel vectors is not fixed at the classical level even when 
one includes a finite J\, since the total exchange field exerted 
at a given site of one sublattice from the four neighboring 
spins of the other sublattice adds up to zero. By contrast, the 
classical energy per site of the 4-spin exchange term goes as 
+KS 4 (cos 2 6> + 2) /6, and thus the plaquette term favors the 
orthogonal state with 6 = ir/2. However up to this point we 
have neglected quantum fluctuations. At sufficiently small x, 
Ji will dominate over K, and will drive quantum fluctuations 
which at the harmonic level are known 20 to favor one of the 
two collinear phases with 9 = or ir (see Fig. 2(b)). Hence 
there must be a critical value of x below which the orthogonal 
state becomes unstable towards the Z 2 collinear phase. 

It is should be noted here that out of the three members of 
Eq. (27), it is the combination JCh — \K. V that is responsible 
for the selection of the orthogonal phase at the classical level, 
since the energy of the third term JC X does not depend on the 
angle 9. 



VIII. NON-LINEAR SPIN- WAVE (NLSW) THEORY IN 
THE EFFECTIVE MODEL 

A LSW expansion around the orthogonal phase in the ef- 
fective model does not capture the instability to the collinear 
phase at small x. This can be seen in the lower panel of Fig. 
4 , where we have also included (by green diamonds) the the 
results from such a calculation (specific details are again rele- 
gated to App. A)) for the correction SS of the 4-fold sites. We 
see that 5S does not show any anomaly down to x = (where 
it approaches the same limiting value with that obtained by 
LSW theory in the full Cairo lattice model). 

A LSW expansion around the collinear phase does not work 
either, since the collinear phase is not the classical minimum 
and the resulting Hamiltonian matrix is not positive definite. 
To stabilize the collinear state, one is then lead to include an- 
harmonic corrections to the theory. This situation is analogous 
to the case of the triangular AFM in a field, where the up-up- 
down state is not the classical minimum but is stabilized by 
the leading 1/S corrections to the linear theory. 37 

Here we have performed a similar Non-Linear Spin- Wave 
(NLSW) expansion around the collinear phase. The quartic 
terms are treated by a standard mean-field decoupling based 
on the following (real) pairing fields which must be finite in 
the collinear phase 

^0 = (at at) = (f+fi) 

Ni = (atfi) = (atfi +x+y ) 

L\ = (difi± x ) = (a>ifi±v) (31) 

L 2 = (ciia i±x ) = (a % a l±y ) = (f t f i±x ) = (fifi± y ) ■ 

Here No is the onsite correction to the spin length, Ni is 
the short-range correlationalong the FM lines of the collinear 
phase, L\ is the short-range correlation along the AFM lines 
of the collinear phase, and finally L 2 is the AFM correlation 
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FIG. 10: The unit cell of the Cairo pentagonal lattice along with 
a clarification of the space group symmetries. The vectors and 
t y denote the primitive translations of the Bravais square lattice. In 
addition to the C4 rotations around the 4-fold coordinated sites we 
also have the four non-symorphic operations (<Ti|T) which stand for 
reflections Oi ii — x, y, d, and d') followed by the non-primitive 
translation t. In the absence of the 3-fold sites, this point group 
reduces to Cav which is the point group of the square lattice. 



between next-nearest neighbors. The mean-field decoupling 
leads to a quadratic theory of the same form with Eq. (Al) and 
with the matrices and Dfe given explicitly in App. A4. 

Given the presence of the two inter-penetrating Neel sub- 
lattices, which is dictated by the dominant J2 coupling, we 
expect that L2 is finite even at x = 0, while both Ni and L\ 
should vanish at x = and become finite as soon as the order- 
by-disorder effect takes place. Our results from the numerical 
self-consistent solution for the above pairing fields are shown 
in the left panel of Fig. 9. As expected, both Ao and L2 are 
finite (and strong) at x = with little dependence on x, while 
Ni and L\ approach zero at small x. 

The right panel of Fig. 9 shows the development of a finite 
gap for the second modes at k = and k = (tt, tt). As we 
show in the inset, these momenta map to the k = (0, tt) and 
k = (tt, tt) mode of the BZ of the J\- J2 model which we get 
if we neglect the K-term at very small x. So the restoration 
of the gap at these k-points is consistent with our expectation 
of having two Goldstone modes in the collinear phase which, 
in the BZ of the J1-J2 model, sit at k = (0, 0) and at k = 
(71", 0) (which is the ordering vector of the state about which 
we performed our semiclassical expansion). The lower panel 
of Fig. 9 shows the dispersion of the two spin-wave branches 
along certain symmetry directions of the BZ. 



IX. EXACT DIAGONALIZATIONS IN THE CAIRO 
LATTICE FOR S = 1/2 

In the remaining part of the article we discuss our exact 
diagonalization results from finite-size clusters with periodic 
boundary conditions and with spin S = 1/2. We have in- 
vestigated both the Cairo lattice model, as well as different 
variations of the effective model in the square lattice. 

The main results from the Cairo clusters can be summa- 



TABLE I: Character table of the little group of k = which is iso- 
morphic to d v . The decomposition of the 5 IR's (first column) 
into modes with well defined angular momenta I are given inside 
the parentheses. 





E 


c 2 


d 


(ff a |r) 


Mr) 










KM 


Mr) 


Ai (I = 0) 


1 


1 


1 


1 


1 


A 2 (I = 0) 


1 


1 


1 


-1 


-1 


Bl {I = tt) 


1 


1 


-1 


1 


-1 


B 2 (I = tt) 


1 


1 


-1 


-1 


1 


£(Z = ±f) 


2 


-2 












rized as follows. At x ~ 2, we find a GS levelcrossing to 
a state with total spin A/6, which is the onset of the 1/3- 
ferrimagnetic phase (Sec. IX B). In the opposite regime of 
small x (Sec. IXC), we find a singlet GS and a low-energy 
structure which proves the presence of two nearly decoupled 
sublattices, in agreement with the effective model. This is 
explicitly demonstrated by comparing the spectra of the 24- 
site Cairo cluster to that of the effective 8-site square cluster 
(Sec. IXC 1). We also examine a number of GS properties 
at small x, such as spin-spin, dimer-dimer, and vector-chiral 
correlations (Sec. IX C 2), as well as the low-energy symmetry 
properties of the spectrum (Sec. IXC 3). These results show 
a strong competition between the collinear and the orthogonal 
state at low x, and highlight the fact that the locking between 
the two sublattice Neel vectors (favored by the large J 2 / J\ in 
the effective model) takes place at large length scales. 

At intermediate x (x ~ 1.2 for the 24-site cluster) we find 
a GS level crossing to a new singlet state, accompanied by 
a whole reorganization of the low-lying excitations and their 
symmetries. This intermediate phase is discussed separately 
in Sec. X. 



A. Finite-size clusters and symmetries of the Cairo lattice 

As we mentioned above, the Cairo pentagonal lattice has 
a unit cell of six sites. It turns out that the largest finite-size 
cluster (with periodic boundary conditions) that has the full 
point group symmetry of the infinite lattice and which is also 
accessible by our computational capabilities has 24-sites (the 
next symmetric cluster has 48 sites which is too large). As we 
are going to see, most of the valuable information presented 
below comes from this 24-site cluster, exactly because it has 
all point group symmetries. We have also investigated clusters 
with 12, 18 and 30 sites which are however either too small or 
lack some of the point group symmetries of the infinite lattice. 

Before we present and analyze our ED results it is also use- 
ful to understand the space group structure of the Cairo lat- 
tice. Its Bravais lattice is the square lattice with primitive 
translations t x and t y along the x- and y-axis respectively 
(see Fig. 10). Apart from primitive translations, the lattice is 
also invariant under the C4 rotations around the 4-fold sites as 
well as under four non-symorphic "glide" operations (<7i|r) 
(i = 1-4), which are reflections followed by the non-primitive 
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FIG. 11: Left panel: Low-energy spectrum in the 24-site Cairo pentagonal cluster as a function of (x — l)/(x + 1) (bottom x-axis), where 
x = J43/J33 (top x-axis). The energy is measured from the respective GS and is scaled with the function f(x) = J33X 2 for x < 1 and 
f(x) = J43 for x > 1. Open symbols stand for S z = states while crosses denote S z > states. The legends specify the different symmetry 
sectors for this cluster. The last portion of the symmetry label accounts for the parity under spin inversion, with "Sze" and "Szo" specifying 
respectively even and odd sectors. The arrow denotes the 1/3-ferrimagnetic state (here with total spin S = 4) which becomes the GS for 
x > 1.96. Right panel: Low-energy spectrum of the 24-site Cairo pentagonal cluster versus the total spin S for x — 0.8. Different symbols 
correspond to different symmetry sectors as shown by the legends. The parentheses next to the legends denote the degeneracy of each sector 
(apart the Zeeman degeneracy). 



translation r (see Fig. 10). Thus the point group Go of the 
Cairo lattice is isomorphic to C 4 „. 

Next, we would like to discuss the Irreducible Represen- 
tations (IR's) of this space group with emphasis on special 
points in the Brilloin Zone (BZ). We begin with the zero mo- 
mentum sector. The little group of k = is the full point 
group and its IR's are taken over from those of C 4 „. These are 
shown in Table I where we also show their decomposition into 
IR's of the C4 subgroup which are labeled by the angular mo- 
menta /. Specifically we have four one-dimensional IR's out 
of which the first two, A\ and A2, are s-wave states (I = 0) 
and the remaining ones, B\ and B2, are d-wave states (I = tt). 
In addition there is a 2-dimensional sector "E" which decom- 
poses into I = ±7r/2 and will be denoted in the following by 
"0.{±tt/2}". 

The little group of k = (7r,7r) is also C± v . However the IR's 
are now different from the k = case, since the representa- 
tion theory for non-symmorphic groups is more involved for 
momenta that sit on the BZ boundary. It turns out 62 that there 
are only two 2-dim IR's for k = (tt, tt). The first combines 
I = and I = tt while the second combines I = ±7r/2. So 
in the following we shall denote these IR's by "(tt, tt).{0, tt}" 
and "(tt, 7r).{±7r/2}" respectively. 

We finally discuss the (tt, 0) point. Here the little group is 
isomorphic to Civ but the IR's are not those of 6*21, since we 
are again dealing with a point on the BZ boundary. Here group 
theory predicts a single 2-dimensional IR 63 whose members 
have angular momenta I = and tt. By including the second 
member (0, tt) of the star of k, we get a single 4-dimensional 
IR which we shall denote here by "(tt, 0)*.{/ = 0, 7r}". 



B. Large-a? regime: The 1/3 ferrimagnetic phase 

Our ED results of the full model show that the 1/3 ferrimag- 
netic state is indeed stabilized at large x. The transition to this 
state can be easily identified by a GS level crossing between 
the lowest 5 = state and the lowest S = N/6 state. For all 
clusters investigated the transition occurs around the classical 
point, x ~ 2. For the 24-site cluster which is the most sym- 
metric cluster, the transition takes place at x ~ 1.96 (see left 
panel of Fig. 1 1 below). 



C. Small-z regime: The presence of two weakly coupled AFM 
sublattices 

1. Structure of spectrum at x — > 

Figure 11 shows the low-energy spectrum of the 24-site 
Cairo lattice as a function of (x — l)/(x + 1). To highlight 
the dominance of the J2 energy scale at small x and to bet- 
ter examine the resulting splitting of the x = GS manifold, 
we have plotted the energy (always measured from the GS) in 
units of x 2 J33 for x < 1, and in units of J43 for x > 1. 

The spectrum shows the following features in the limit of 
x = 0. The first excitations above the singlet GS are two 
triplets with excitation energy x 2 J33. At energy 2x 2 Jzz we 
find all together nine degenerate states, among which we have 
three singlets, five triplets and 1 quintet. 

To understand this structure we turn to the effective model 
description of the 24-site Cairo pentagonal cluster (see left 
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FIG. 12: Left panel: The effective J1-J2-K model description of the 24-site Cairo pentagonal cluster. There are two inter-penetrating square 
plaquettes which consist of the sites (1, 0, 3, 5) and (2, 4, 7, 6). The Ji and J2 interactions are denoted by thin dotted and thick solid lines 
respectively, while the orientation of the A"-term in each Ji-plaquette is denoted by a thick segment (which coincides with the underlying 
J;!3-dimer of the Cairo lattice) at the middle. At small x the two square plaquettes are decoupled from each other and only the J2 coupling 
survives (here J = 2 J2). The total energy spectrum (last panel) results from a direct sum of the energy spectra of the two plaquettes (second 
and third panels). In addition to the spin quantum numbers we also specify the angular momentum I with respect to C4 rotations around the 
site numbered 5. All degenerate levels are split by hand for clarity. 



panel of Fig. 12). This cluster comprises 8 four-fold coordi- 
nated sites which form two inter-penetrating square plaquettes 
A (sites 1, 0, 3, and 5) and B (sites 2, 4, 7, and 6). This sim- 
plifying aspect allows to work-out explicitly the low-energy 
structure of the effective model for small x. Let us label the 
states of plaquette A by |So5'S , 13'S'a), where S05 = So + S5, 
S13 = Si + S3, and Sa = S05 + S13. Similarly the states of 
plaquette B are labeled as \S27S4qSb)- To lowest order in x 
only the J2 terms survive which read 

7^2 = 2 J2 ^S4g • S27 + S05 • Si3^ 

= J-2 (s\ + S% — S 2 6 — S27 — Sq 5 — S 2 3 ^) , (32) 

where the factor of two in the first line accounts for the fact 
that each J2 interaction appears twice in our cluster due to 
the periodic boundary conditions. Thus the two plaquettes 
A and B are decoupled from each other, and the full spec- 
trum E = Ea + Eb (rightmost panel of Fig. 12) can be 
obtained by adding the two single-plaquette spectra (second 
and third panels of Fig. 12) by a standard addition of angular 
momenta S = S 4 + Sg. It is then straightforward to show 
that the global GS is the direct product of the two 1 1,1,0) 
singlets of each plaquette and that it is s-wave (I = 0). The 
lowest excited states are the triplets |1, 1, 1) a ® |1, 1, Q)b and 
|1, 1, 0)a ® |1, 1, l)s, which have energy J = 2J 2 above the 
GS and angular momentum I = and it respectively. At en- 
ergy 2 J, we get 3 singlets, 5 triplets, and 1 quintet. 

So the multiplicities and the symmetry properties of the 
spectrum at small x match exactly the ones found by ED in 
the original 24-site Cairo cluster (Fig. 11). This confirms the 
main picture from the effective model of having two nearly 
decoupled inter-penetrating square AFM's. 



2. GS correlations 

GS spin correlations. — Figure 13 shows the GS spin-spin 
correlations for the 24-site cluster. The left panel shows the 



correlation profiles at a representative value x = 0.5, while 
the right panel shows the dependence as a function of x. There 
are two important features in these figures. First, there is a 
strong AFM correlation between the reference site and the 4- 
fold sites that belong to the same subsystem. This confirms 
the strong AFM J2 coupling within each sublattice. 

The second feature, which remains true in the entire range 
up to x = 1, is that there are almost no correlations between 
the reference site and the 4-fold spins belonging to a differ- 
ent sublattice. This feature is consistent with the orthogonal 
phase scenario where the two Neel vectors are perpendicular 
to each other (9 = n/2). It would also be consistent with the 
collinear phase scenario where the correlations actually vanish 
for finite-size clusters because the GS is an equal superposi- 
tion of the two genuinely different collinear phases (6 = and 

7T). 

Hence the spin-spin correlation data cannot establish the 
critical value of x where the transition between the orthogo- 
nal and the collinear phase takes place. For this, one would 
need much larger cluster sizes since, as we know from pre- 
vious studies on the J1-J2 model, 27 ' 28 there is a very large 
length scale associated with the locking between the two Neel 
vectors. 

GS dimer and vector-chiral correlations. — We now turn to 
some other correlations which may also be used as diagnos- 
tic tools for the collinear and the orthogonal phases. In the 
collinear phase the spins order FM in one direction and AFM 
in the perpendicular direction, so the dimer-dimer correlations 
along the two directions should reflect this physics. 28 On the 
other hand, in the orthogonal phase the spins lie in one plane 
so the orthogonal phase must show a staggered signal in the 
vector-chiral correlations. 31,32 Figure 14 shows these two type 
of correlations for the same GS at a representative value of 
x = 0.1. The results show that both types of correlations are 
present in the same GS and with the expected profile. As in 
the case of the spin-spin correlations, this shows that the two 
phases compete with each other but the 24-site cluster is too 
small to discriminate between the two. 
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FIG. 13: Left: Ground state spin-spin correlation profiles (Se • Si) at 
x = 0.5 for the 24-site Cairo pentagonal cluster. Filled (blue) circles 
denote positive correlations while open (red) circles denote negative 
correlations. Right: Correlations between the reference site (6) and 
the set of all symmetry inequivalent sites of the cluster as a function 
of x, up to x = 1. Note that the limiting correlation values at x = 
are consistent with the exactly known GS wavefunction of the 8-site 
effective cluster (see text). 



3. Excitations: Low-energy Towers of States 

We next turn to the low energy excitations and check 
whether we see any signatures of the collinear and the or- 
thogonal states. This is done by looking at the symmetry 
structure of the so-called Anderson towers of states. 38 It is 
by now established, following the seminal work by Bernu 
et al. 39 ' 40 and Lecheminant et al. 41 ' 42 , that a given magnetic 
phase in the thermodynamic limit shows up in finite-size spec- 
tra through the clear formation of a tower of states which scale 
as S(S + 1)/N and is well separated from higher excitations. 
A wavepacket out of this infinite tower would be stationary in 
the thermodynamic limit and would correspond to the given 
classical state. Not surprisingly then, the multiplicities and 
symmetry properties of this set of states are intimately con- 
nected to the symmetries that are broken in the classical phase 
and can actually be derived by group theory alone. 39 ^ 4 

Now, the collinear and the orthogonal phase break the full 
symmetry group of the Hamiltonian in a different way, so the 
structure of the corresponding tower of states should be very 
different from each other. In App. F we derive the symmetry 
content of the two towers using group theory. 

The final predictions for the lowest total spin S sectors of 
the towers are given in Tables IV and V. Specifically, the 
collinear phase comes with 2 states per total spin S, while 
the orthogonal phase should show 25* + 1 states (not related 
to the Zeeman degeneracy) at a given spin sector S. The sym- 
metry properties of these states with respect to C4 rotations as 
well as the four non-symorphic operations (er^T) are shown 
in Tables IV and V. 

These towers should be now compared to the low-E excita- 
tions of the 24-site cluster which are shown in the right panel 
of Fig. 1 1 as a function of the total spin S, at a representative 
value x = 0.8. In the singlet sector we find the "0.A1" ground 
state, a "0.B1" state nearby and another "0.B2" state slightly 
higher in energy. The first one belongs to the orthogonal tower 
but we may also think of the pair "0.A1" and "0.B2" as parts 
of the collinear tower. The "intruder" state "0.B1" does not 
belong to any of the towers and indeed at small enough x this 




bond nematicity vector chirality 

FIG. 14: Ground state bond nematic ((Si ■ S2)(Sk ■ Si)) and vector- 
chiral ((Si x S2) z (Sfc x S;) z ) correlations at x = 0.1 for the 24-site 
Cairo pentagonal cluster. The thick bond (1,2) is the reference bond. 
Solid (blue) lines denote positive correlations and dashed (red) lines 
denote negative correlations. The width of each line is proportional 
to corresponding expectation value. 



state is slightly higher than the "0.B2" state. At x = 0.8 this 
state has very low energy because it is about to become the 
GS above x ~ 1.2 (see Sec. X below). 

Likewise, in the 5 = 1 sector, we find three triplets very 
close in energy. Two of these belong to the two-fold IR 
"(■7T, 7r).{0, 7r}" expected for the collinear tower, but we can 
also include the third state "0.B1" to complete the three states 
expected for the orthogonal tower. A similar situation occurs 
in the S = 2 sector. Altogether, we find that the low-E exci- 
tations contain both towers of states which conforms with the 
previous picture from the GS correlations, that the collinear 
and orthogonal phases compete with each other but the 24- 
site cluster is too small to discriminate between the two. 



X. THE INTERMEDIATE NON-MAGNETIC PHASE. 

A. Identifying the relevant effective term which drives the 
intermediate phase 

We now turn to intermediate values of x. Looking back at 
the low-E spectra shown in the left panel of Fig. 1 1 we see 
that there is a GS level crossing in the singlet sector around 
x ~ 1.2. Similar GS level crossings are also found in the other 
Cairo clusters that we studied. For the 24-site cluster the new 
singlet GS belongs to the "0.B1" sector which belongs to nei- 
ther the collinear nor the orthogonal towers of states (Tables 
IV and V), which is the first strong evidence that the system 
enters a new phase. 

Another evidence is the reorganization of the low-E excita- 
tion spectrum. In particular as we see in left panel of Fig. 15, 
the energies of the three lowest spin sectors form a concave 
envelope, which suggests a tower of states formed by even 
spin sectors only. Indeed the magnetization (not shown here) 
grows in 5S Z = 2 steps until we reach the 1/3-ferrimagnetic 
phase. Provided that these SS = 2 excitations form a tower 
that collapses in the thermodynamic limit, the resulting state 
would not break time-reversal invariance, as it happens e.g. in 
a spin-nematic state. 

It turns out that we can actually learn more by tracing the in- 
termediate "0.B1" state back to its original place in the x = 
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FIG. 15: Low-E spectra as a function of total spin 5 for three different cases. The left is the spectra of the 24-site cluster for x = 1.667, i.e. 
deep inside the intermediate phase. The middle panel shows the spectrum of the 8-site unconstrained J1-J2-K model in the large K/Ji t 2 
regime. The third panel shows the spectrum of the 8-site effective model but this time we have kept only the K x term among the three plaquette 
terms of Eq. (26). The pairs of numbers inside the parentheses indicate the total spins Sa and Sb of the two sublattices which are separately 
conserved in this SU(2)xSU(2) model. In all three cases, the ovals indicate the sequence of GS's visited by the system in a field, while the 
different symbols follow the same convention as in the right panel of Fig. 11. 



spectrum. First of all, Fig. 11 tells us that the "0.B1" state 
is a member of the x = GS manifold where J 33 -dimers 
form singlets, in contrast to the 1/3-ferrimagnetic state which 
is clearly not. So the intermediate phase is more closely re- 
lated to the physics of the x = limit. Secondly, we know 
explicitly the low spin states of the 24-site cluster at x = 
and so we may find out exactly which excitation at x = 
evolves into the intermediate "0.B1" phase. We find that it 
is the following combination among two of the singlets with 
energy 2 J above the GS (see right panel of Fig. 12): 

I*) = -^(|1,1,0)a® 10,0,0)^ + |0,0,0)a® |1,1,0) b ) . 

(33) 

This state has an important property (not shown explicitly 
here): It minimizes simultaneously all plaquette terms of the 
type K x of Eq. (27), with eigenvalue -3/16. This provides ev- 
idence that the instability mechanism that triggers the tran- 
sition to the intermediate phase is related to the four-spin 
exchange term JC X . At the level of the effective model of 
Sec. VII, this suggests that the K x terms dominate in much 
higher orders of perturbation theory in x, but we are not able 
to check this explicitly. 

We can still however make progress using only the three 
couplings that we know so far from the fourth-order theory. 
The idea is to compare the low-E spectrum of the intermediate 
phase of the Cairo model with that of the unconstrained J\- 
J2-K model as we visit different regions in the J1/J2-K/J2 
plane. We have performed ED in the effective lattice model 
(which has 2 sites per unit cell) using clusters with 8, 10, 16, 
18, 20, 26, and 32 sites. The fully symmetric clusters with 8, 
16 and 32 sites give the most clear and systematic evidence 
so we shall only discuss these clusters here. We begin with 



the 8-site effective cluster since it can be directly compared to 
the 24-site Cairo cluster (the 4-fold sites of the latter form the 
8-site effective cluster). By an inspection of the low-E spectra 
in various regions of the J\j J2-K/J2 plane, we have located 
the region of large Kj J1.2 as the one with very similar low-E 
spectral features with that in the intermediate phase. This is 
demonstrated by the first two panels of Fig. 15. 

Given now the special role of K x regarding the above state 
|$) (see Eq. (33)), we next check what happens if among the 
three terms included in K (see Eq. (26)) we keep only the JC X 
term. The corresponding 8-site spectrum is shown in the third 
panel of Fig. 15 and demonstrates that the spectrum retains 
the same features at low energies whether we keep only JC X or 
not. This shows that K x is indeed the most relevant plaquette 
term for the intermediate phase. 

B. The K x -model 

For the above reasons we shall restrict ourselves to the JC X - 
model in the following. The Hamiltonian reads 

(S« • S 7 ) (S^ • S d -) (34) 

plaquettes 

where the sum is over all plaquettes of the square lattice, and 
(a, /?, 7, 6) label the spins around a plaquette clockwise. This 
model can be thought of as a 2D generalization of the well- 
studied 45 ^ 9 ID spin-orbital model. 

We first discuss the classical limit of this model. It is easy 
to show (and we have also checked it by Classical Monte 
Carlo simulations) that the energy can be minimized by any 
collinear configuration with "3up-ldown" or "3down-lup" 
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FIG. 16: (Color online) First two panels: Low-E spectra as a function of total S z for the 16- and 32-site clusters using only the effective K x 
plaquette term. Here all energies are measured for the GS energy Eo(S z ). The ovals indicate the sequence of GS's visited by the system in 
a field while the pairs of numbers inside the parentheses indicate the total spins Sa and Sb of the two sublattices in these states. Note that 
for the 32-site case, the 1/2-magnetization states with (Sa, Sb) = (0, 8) and (8, 0) are not the GS's. Right panel: Magnetization process in a 
field for the clusters with 8, 16 and 32 sites. 



spins in every plaquette. In particular, starting from one such 
GS we can generate others by flipping all the spins along any 
horizontal or vertical line of the square lattice. This leads to 
a sub-extensive number of GS's (with variant magnetizations) 
which are not related by global spin rotations. So the classical 
limit is highly frustrated and thus we anticipate rich physics 
from thermal (not studied here) or quantum fluctuations. 

It is natural that the above classical GS manifold appears 
also if we use Ising spins in Eq. (34). Actually, as found by 
Rojas et al., m the same GS's arise also in the Ising version of 
Eq. (1). 

Turning now to the quantum limit, let us denote by A and 
B the two interpenetrating square sublattices of the full square 
lattice (in Fig. 18 of App. E, the A and B sublattices consist 
of the (a, 7)- and (/3, <5)-bonds respectively). This model has 
a C4 rotation symmetry around the center of the plaquettes, 
contrary to the full J1-J2-K model. More importantly, the 
model has an enlarged SU(2)xSU(2) symmetry since we can 
make independent spin rotations in the two sublattices A and 
B without changing the energy. So the total spins Sa and Sb 
of the two sublattices are good quantum numbers. This also 
means that within a fixed (Sa, Sb) manifold the energy is 
independent of the total spin S. In addition, the model has an 
Z2-invariance under interchanging the two sublattices A and 
B, so the manifolds (Sa, Sb) and (Sb, Sa), with Sa 7^ Sb, 
must be degenerate. 

Let us now try to establish some systematic spectral fea- 
tures in the JC X -model by looking at the larger 16- and 32-site 
symmetric clusters. Their spectra are shown in the first two 
panels of Fig. 16 as a function of total S z . All energies are 
now measured from the GS energy Eq(S z ) so that we better 
retrieve the details of the spectra. The states that are visited by 
the system in the presence of a magnetic field are highlighted 



by ovals (tower of states in the following), while the pairs of 
numbers inside the parentheses show the total spins (Sa, Sb) 
in these states. 

All towers of states, including the one for 8-sites (third 
panel of Fig. 15) exhibit common features. Apart from the 
singlet GS with (Sa,Sb) = (0,0), and the fully polarized 
state (N/4, N/4), all other states show a two-fold degener- 
acy which is related to the Z 2 -symmetry mentioned above. 
More importantly, there is a very clear difference between the 
states below and above 1/2-magnetization. The tower below 
1/2 comprises only states with even S, while above 1/2 we 
have SS — 1 steps. This is also demonstrated in the third 
panel of Fig. 16 which shows the magnetization process in a 
field. Altogether these results reveal two different states above 
and below 1/2-magnetization. 

The physics above 1/2-magnetization can be easily under- 
stood by noting that the corresponding states have either Sa 
or Sb equal to N/4. Thus one of the two subsystems is fully 
polarized above the 1/2-magnetization. It is straightforward 
to show that in this case the /C^-model reduces to the two- 
spin exchange model in the other sublattice with a nearest- 
neighbor coupling equal to K x /4. We have checked that the 
tower of states above the 1/2-magnetization match exactly 
(both in the symmetries and the actual energies) the corre- 
sponding tower of states of the square lattice AFM with half 
the number of sites and a nearest-neighbor coupling equal to 
K x /4. 

The physics below 1/2-magnetization is much more inter- 
esting. We first note that the tower below 1/2 contains states 
where both Sa and Sb are even, and that one of the two is 
always zero. This suggests that the SU(2) symmetry is bro- 
ken down to U(l) only in one of the two sublattices. Namely 
that we have some kind of spin nematic state in one sublat- 
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FIG. 17: (Color online) GS nematic correlations ({Sf Sj){S^ S^)) 
for the /Cz -model in the 32-site cluster. The cluster is denoted by the 
(green) rectangle. The think (black) line is the reference bond (y). 
The solid (blue) lines denote the bonds (fcZ) with positive correla- 
tions while the dashed (red) lines denote negative correlations. The 
width of the bonds is proportional to the actual strength of the corre- 
lations, the largest of which is equal to 0.0228. Note that there are 
no correlations between the reference bond and bonds belonging to 
a different sublattice. 



tice and a spin liquid state in the other sublattice. In addition, 
the angular momentum alternates between I = (s-wave) and 
I = 7r (d-wave) as we go down in S z starting from the 1/2- 
magnetization state which is s-wave. Hence the zero-field GS 
is s-wave for 16 and 32 sites but it is d-wave for 8 sites, which 
conforms with the d-wave property of the intermediate "0.B 1" 
state in the 24-site Cairo cluster. 

The alternation between s-wave and d-wave GS symmetry 
at low-magnetizations is exactly what happens in the model 
studied by Shannon et al. in Ref. [29], and suggests that 
the intermediate spin-nematic phase has a d-wave symme- 
try. We have checked this numerically by a calculation of the 
GS bond-nematic correlations of the type ((Sf Sj)(S^ <5j~)) 
where and (k, I) denote two different bonds of the lat- 
tice. The results are shown in Fig. 17 and confirm the presence 
of a d-wave signal in one sublattice only. The correlations be- 
tween bonds of different sublattices vanish exactly because 
the GS belongs to the sector (S A , S B ) = (0, 0). 

Let us now try to identify the low-E degrees of freedom that 
could give rise to such a nematic state. Since we are dealing 
with spins 1/2, we need some composite degrees of freedom 
that can sustain a finite quadrupolar order parameter. For the 
spin-nematic phase of Ref. [29], these might be the triplets 
formed on nearest-neighbor bonds, 29 or the plaquette 5 = 2 
states, as discussed by Ueda and Totsuka. 50 In the present case 
we can also identify some triplet degrees of freedom by con- 
sidering the effect of JC X on a single plaquette (a, /3, 7, 8), It 
is easy to see that there are two possible ways to minimize the 
energy, either by forming a singlet on the (a, 7)-bond and a 



triplet on the (/3, <5)-bond or vice versa. Both states have en- 
ergy —3/16 and have total spin S = 1, so they can indeed 
sustain a finite quadrupolar order parameter. 

Now, there are two factors that impose nontrivial con- 
straints into the way these triplet degrees of freedom inter- 
act with each other in the full lattice. The first is that, by 
their nature, the two triplet GSs impose that a triplet in one 
diagonal bond must come with a singlet on the other diago- 
nal bond. The second complication is that any given site par- 
ticipates into four plaquette interactions. One approach that 
deals with the first problem (but not with the second) is to 
perform a strong-coupling expansion around the limit where 
only 1/4 of the plaquettes have a finite K x and are decoupled 
from each other (see Fig. 18). 50 ' 51 The details of this approach 
are provided in App. E and follow closely in spirit the treat- 
ment by Lecheminant and Totsuka 51 of a very similar situation 
in a two-leg ladder system where a different pair of plaque- 
tte triplets emerges at low energies. This approach delivers 
an effective model which is similar to the well-known Kugel- 
Khomskii Hamiltonian 52 for orbital degenerate systems. Here 
the role of the orbitals is played by the two types of triplets. 

A natural variational treatment of this effective model de- 
livers a spin-nematic GS of a novel kind where both types of 
triplets are entangled (see App. E). Despite this success, there 
are two problems with this wavefunction. The first is that 
both sublattices of the square lattice participate in this state 
(see App. F 4 for the symmetry decomposition of this state), 
while the numerical tower of states suggest that the symmetry 
breaking occurs only in one sublattice. The second problem is 
that the variational state breaks translational invariance which 
is also in disagreement with the numerical tower of states. In 
particular, the d-wave character of the correlations is built-in 
from the outset since, by construction, the "plaquettized" lat- 
tice distinguishes between the two different directions in each 
sublattice (see App. E). 64 



XI. SUMMARY 

In this article we have presented an extensive analytical and 
numerical study of the Heisenberg model on the Cairo pentag- 
onal lattice. We have shown that by tuning the spin S and the 
ratio x of the two inequivalent exchange couplings of the lat- 
tice, we can drive the system through a number of competing 
phases, some of which have a strong quantum-mechanical ori- 
gin. We have demonstrated that the rich physics of this model 
arises, to a large extent, from the presence of two inequivalent 
sites in the Cairo lattice. To highlight this central feature and 
to better understand the nature of the various phases we have 
followed a series of different approaches which build quantum 
fluctuations in a gradual way starting from the classical limit. 

The classical phase diagram contains three phases, the or- 
thogonal phase for x < \/2, the 1/3-ferrimagnetic phase for 
x > 2, and a mixed non-coplanar phase for V2 < x < 2 
which combines the orthogonal and the ferrimagnetic. The 
orthogonal state is the configuration found experimentally for 
the Bi 2 Fe 4 Og compound, 13 which suggests that this phase is 
quite robust in a more extended parameter space compared to 
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our symmetric version of the model. 

In the quantum 5 = 1/2 limit, we have found an under- 
lying order-by-disorder mechanism which is active at small x 
and which stabilizes a collinear magnetic configuration. The 
1/3-ferrimagnet phase seems to survive quantum fluctuations 
down to 5 = 1/2 in the large-.T regime. On the other hand, 
the mixed phase becomes unstable for low enough 5 and dis- 
appears completely from the phase diagram for 5 = 1/2. 

Our ED results provide strong evidence of a new intermedi- 
ate phase with spin-nematic correlations. We argued that this 
phase is stabilized by an effective 4-spin plaquette interaction 
that first appears in fourth-order perturbation theory in x. This 
term has an enlarged SU(2)xSU(2) symmetry structure sim- 
ilar to spin-orbital models, and favors two species of triplet 
GSs on a single plaquette. The symmetry structure of the low- 
E spectra in the corresponfing plaquette model suggests that 
these triplet degrees of freedom order in a non-trivial fashion, 
whereby one of the two square sublattices forms a transla- 
tionaly invariant spin-nematic structure with d-wave symme- 
try (similar to the one found by Shannon et al. 29 ), while no 
symmetry-breaking seems to occur in the second sublattice. 

One of the novel aspects of the intermediate phase is its re- 
sponse under a magnetic field. The sublattice where the spin- 
nematic order takes place responds by developing a dipolar 
moment until full saturation. At this point, which corresponds 
to 1/2 moment of the full system, the low-energy physics of 
the second sublattice effectively reduces to that of the square 
lattice Neel AFM. 

A simple physical picture for the intermediate phase is still 
lacking. We have discussed e.g. a strong-plaquette expansion 
that builds upon the low-energy triplet degrees of freedom of 
a single plaquette. A variational treatment of the resulting 
model does deliver a spin-nematic phase, but this seems to en- 
tangle both sublattices and in addition it has symmetry proper- 
ties that are not compatible with our low-energy spectra from 
ED. It is our hope that further studies will shed more light in 
this direction. 

At a more general level, this study demonstrates that pen- 
tagonal lattice models can feature a wealth of strong corre- 
lation phenomena. Compared to the well-explored triangle- 
based models (triangular, kagome, etc.), here we have a qual- 
itatively different degree of frustration and complexity since 
the low-energy physics is not related to the local single- 
pentagon physics in any obvious way. Hence pentagonal lat- 
tice models are interesting in their own right and provide a 
distinct platform for exploring and realizing novel phases of 
matter in frustrated magnetism. 
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Appendix A: Spin Wave Theory (LSWT) 

1. General procedure 

Here we outline the general procedure for diagonalizing the 
quadratic bosonic Hamiltonians obtained from LSWT around 
the three classical ground states. The Cairo pentagonal lat- 
tice has 6 sites (two 4-fold and four 3 -fold coordinated) per 
unit cell. So we introduce six bosonic operators denoted by 
a,i,bi,Ci,di,ei, fi to describe the harrmonic fluctuations on 
each site of the i-th unit cell. Quite generally, up to quadratic 
order the semiclassical expansion gives 



n = E c 



5E 1 



Hi. • Ai 



(Al) 



where SEi is a constant (see below), 



A t = { a k' b k> c t> d k> e k>fk> a -ki b -k,c-k,d_ k ,e_ k ,f- k ) 

(A2) 

and Hfc = I _ + _ T I , where both Cfe and are 6 x 6 

V U k C -fc / 

matrices. To diagonalize the Hamiltonian we search for a 
new set of bosonic operators A k given by the generalized Bo- 
goliubov transformation Afe = Vfe • Afe, such that the ma- 
trix V^HfcVfc = fl k becomes diagonal. The transformation 
must also preserve the bosonic commutation relations, which 
can be expressed compactly as g = g = V& • g • VjJ", where 
g is the "commutator" matrix 



g = A fc -A+- (At) T -Al 



1 6 

o -i 6 



(A3) 



and lg stands for the 6x6 identity matrix. The above two 
conditions give 



( g H fc ) • v fc = v k ■ (gO fc ) = v fe • n' k 



(A4) 



which is an eigenvalue equation in matrix form (the columns 
of Vfe contain the eigenvectors of gHfe). 

One can further show 34 that if Hfc is semi-definite positive, 

then fife = [ ^ ^ ) , where u> k is a diagonal matrix with 
\ U u) k J 

non-negative entries LU\ k -ujQ k . This in turn leads to 

H = E cl +SE 1 +SE 2 

+ S^2 (uikata k +...+ U6kfkfk) , (A5) 

k 

with 8E-2 = ■§ ^2 k (uJik +• • •+ ^6fe), which represents the 
total zero-point energy from all harmonic oscillators in the 
theory. The total quadratic correction to the GS energy is then 
given by SE = 5E\ + 5E 2 - 

We now turn to the quadratic correction to the local spin 
lengths. Let us consider the spin S a operator inside the unit 
celH = (0,0,0). We have 



SI — S - a+ =0 a,i =a - S - — — ^ a t a q > 

K,q 



(A6) 
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where N uc = N/6 is the number of unit cells. Using the 
transformation = Vfc • Ak we get 

S* = S - J- Y,Y, Vk{l,n)*V q {l,m) A+(n)A k (m) . 

UC k,q nm 

(A7) 

In the vacuum GS the only non-vanishing expectation values 
are of the type (akat) = 1. Thus from the above sum we may 
keep only terms with n = m = 7 — 12 and k = q. Namely 



s*HS--^^f> fe (i,n)| 2 , 



k n=l 

and similarly for the 3-fold sites: 

12 



(Sg) = S--^£|V fc (2,n)f 



(A8) 



(A9) 



2. LSWT in the Cairo pentagonal lattice 



In the following we provide the explicit expressions of the 
6x6 matrices and defined above. We shall make use 
of the following definitions: k xy = k x + k y ,q± = |(1±^=), 

z ± = MiE^M ancl x± = | ± 1. 



In the orthogonal phase E c i = Ne ort ho , SEi = Ne or tho/S, 
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For the 1/3-ferrimagnetic phase we find E c i = Nef en i, 8E\ = N£f en i / S, and 



(A10) 
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Finally, in the mixed phase E c [ = Ne^ed, SEi = Ne m i xei i/S, and 

„2 
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(All) 



(A 12) 



(A 13) 



r 



3. LSW Theory around the orthogonal phase in the effective 
model 



the dimension of the corresponding matrices Ck and is 
equal to 2. For the theory around the orthogonal phase these 



The unit cell of the effective model has two sites whose 
corresponding bosonic operators are labeled by and fi. So 
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matrices are given by 



Cfc = 



/ 4( J 2 - K) Jip(fc) 

V Jip(-k) i(j 2 -k) 



and 



2( J 2 - K)r{k) Jip(fc) 
Jip(-fe) 2( J 2 - K)r(k) 



(A14) 



(A15) 



where AT = KS 2 , r(k) = cosk x + cosk y , and p(k) 



1 



D i(fe x +fc H )~l 



/2. Finally, the classical GS 



energy per site is e' orth J S 2 
-2 J 2 + 2A. 



-2J 2 + A', while ^ 



4. Non-linear Spin Wave Theory around the collinear phase of 
the effective model 



The explicit forms of the 2x2 matrices Cfc and Dfc that 
enter the mean-field decoupled quadratic theory of Sec. VIII 
are the following: 



4J 2 +4A(A-2)+£i/S 
Ji + If (1 - A) -K2/S) (1 + e" 



2(K -J 2 + &/S) cos k x + 2( A -J 2 + cos k y 



( Ji + K(l - A) + ea/s) (1 + e^y ) \ 
4J 2 + AK(\ - 2) + J ' 

- Ji + A(l - A) + fc/s) (e ife * + e lfc ») 



' x + e l,=1 



) 2(K -J2+ &/S) cos + 2{K -J 2 + fr/S) cos fcj, 



6 



Ji + A(l - A) + &/S) (e" 
where we have introduced the parameters 

47V (3(2 - A) A - J 2 \ - 2N 1 {2,k{\ - A) + J x ) - 2L 1 (zK{l - A) — Ji ) + 4L 2 (J 2 - 3A) 
N x ((9 - 5A)A + Ji) - A (3A(1 - A) + Ji) - 4LiA - 2(1 - \)KL 2 
Li((9 - 5A)A - a) - ^(3^(1 - A) - a) - 47ViA - 2(1 - X)KL 2 
L 2 (k(5 - 2A) - J 2 ) + N (A - 3A) + 2AALi - 2AAi 
L 2 (A(5 - 2A) - a) + A (A - 3A) + 2AAAi - 2ALi . 

I 



(A16) 



(A17) 



(A18) 



Appendix B: QM problem of a single AFM dimer in a staggered 

field 



The minimization of the variational ansatz described in 
Sec. V showed that the 4-fold spins remain coplanar and or- 
thogonal to each other up to x = 2. In this configuration the 
two exchange fields that are exerted on the two sites of each 
J33-bond are antiparallel and have magnitude h s = \/2Sx 
(where we have used a length S for the classical 4-fold spins). 
One then realizes that this problem is equivalent to that of an 
AFM dimer in the presence of a staggered field. A similar sit- 
uation appears for the experimental compound Cu 2 Cd(B03) 2 
inRef. [53]. 

The Hamiltonian of an AFM dimer in a staggered field h s 
is given by 



T~L — J Si • S 2 — /i s (S*j — S*2 



(Bl) 



Importantly the staggered field does not commute with the ex- 
change interaction, and therefore it can polarize the system 
immediately. This is in contrast to the case where we have a 
uniform field, where one must exceed the singlet-triplet gap J 
to polarize the system. The triplets \t\) = | "ft) and = 
1 14-) are eigenstates of T~L with energy J / 4, while the singlet 
I a) = U) ^ } , and the triplet \t ) = are admixed as 



follows: H\t ) = { \t )-h s \s},<mdU\s) = -^-\s)-h s \t ). 
A straightforward diagonalization in this manifold gives the 
following eigenstates and eigenvalues 



IV'i) 

|V> 2 ) 



u|«)+v|tn>, ei = - J/4 - ^/j 2 /A + h 2 (B2) 
«|*>-u|to>, e 2 = -J/4 4-/774+7^(^3) 



with u = cos 9, v = sm9, and tan(20) = 2h s /J. As ex- 
pected, ->• |s) for h s -> 0, while for h s > J, -> | tl>- 
The GS expectation values of the local polarizations and the 
exchange energy are given by 



(Vi|Si-S 2 |Vi) 



3 , 



(B4) 
(B5) 



Appendix C: Quadratic fluctuations around the variational GS 

Here we provide the explicit form of the 8x8 matrices Cfc 
and Dfc that appear in Eq. (17). The various constants that 
appear below are defined as follows: q = xySM/4, y± = 
q{v±u),£± =l±l/v / 2. 
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Appendix D: Strong coupling expansion at small x = J43 / J33 



Here we provide some details of the degenerate perturba- 
tion theory around x = discussed in Sec. VII. For x = 0, 
the nearest-neighbor 3-fold sites form isolated singlets and the 
4-fold sites are free to point up or down. For finite x, the 4- 
fold sites begin to interact with each other through the virtual 
fluctuations of the J33-bonds out of the singlet manifold. The 
resulting effective Hamiltonian acts on the 4-fold sites only 
and can be obtained by degenerate perturbation theory around 
x = 0. This expansion can be organized in a linked cluster 
way as for non-degenerate perturbation theory. 54 The effec- 
tive Hamiltonian is a sum of contributions from a sequence of 
finite clusters of sites. For each given order of perturbation 
theory only clusters up to a given size may contribute. 

More explicitly, an 77-th order process that lives on a given 



cluster g can be obtained by 

rH hr n _i=n— 1 

H$= c {r} VVTZ {ri) V ■■■VTZ {r - l) VV . 

ri,...,r n _i>0 

(Dl) 

Here V is the projection in the unperturbed manifold, Q = 
1—V, TZq = (Eq — QHoQ)~ is the resolvent operator, Eq is 
the unperturbed energy, TZ W = -V, while TZ^- 1 ^ = (TZ ) r . 
Finally, the coefficients C{ r } can be found tabulated (up to 
sixth order) in Ref. [55]. 

We have performed the above expansion up to fourth order 
in x using all relevant clusters. In order to avoid having to 
subtract processes that live on the subclusters of g, 54 we en- 
force that we only keep processes that invoke all elements of 
the cluster g. The first column of Table II shows all clusters 
that generate a finite interaction between the 4-fold sites (writ- 
ten explicitly in columns 2-4). By including all possible ways 
that we can embed each cluster on the lattice one obtains the 
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TABLE II: Local effective terms that are generated up to fourth order in the strong coupling expansion in x — J43 / J33. Only the nine clusters 
that give finite terms are shown (cf. first column). Dashed (solid) lines denote J43 (J33) couplings. The last three columns show the total 
contribution of each cluster to Ji, J2 and the amplitude K of the plaquette term of Eq. (27). 
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corresponding contribution to the effective couplings J\, J2 
and K (columns 5-7). 

Polarization on the 3-fold sites: — The above effective de- 
scription lives in a projected Hilbert space where pairs of 
n.n. 3-fold sites pair-up forming exact singlet wavefunctions. 
However the true GS of the problem has also a non-vanishing 
component on the orthogonal manifold. This component must 
be taken into account if we want to find e.g. the polarization 
on the 3-fold sites. Specifically, as soon as the 4-fold sites 
order magnetically they will exert a finite exchange field on 
the 3-fold sites. For example, in the orthogonal phase the to- 
tal exchange field will be staggered in all dimers and thus it 
will immediately admix a triplet component into the singlet 
GS and give rise to a finite (staggered) polarization. This is 
in contrast to the uniform field case where one must exceed a 
critical value (the singlet-triplet gap J) in order to polarize an 
AFM dimer. 



The component of the GS wavefunction out of the singlet 
manifold Q\^f) can be expressed as 

Q\V) =KVV\V}, (D2) 

where K = (E - QUO.)' 1 is the full resolvent. Thus 

\^) = (1+TZV)\V^>) ^ ^(l+TZV)^"), ( D3 ) 

where |^°) is the normalized GS of the effective model 
(times the product of J33 -singlets), and zq = 1/(1 + 
(xEf 1 T^T?. 2 M/ ) ) plays the role of a "wavefunction renormal- 
ization factor" (specifically 1 — zq measures the degree of ad- 
mixture inside the GS from states outside the unperturbed GS 
manifold). 

Now, the GS expectation value of the magnetization say at 
site 1 of a given dimer is then given by 

(S?) = 2 z Re (^I^WI* ). (D4) 
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To lowest order in V we may replace zq ~ 1 and 72. ~ TZq, 
which leads to the standard expression from linear response 
theory. In particular, the result coincides with the linear or- 
der contribution found for the problem of an AFM dimer in a 
staggered field studied in App. B. 

We should remark here that in the collinear phase the total 
exchange field vanishes in half of the J33-dimers, see e.g. the 
site labeled by Cj in Fig. 2(b). So the polarization in these 
dimers must grow quadratically and not linearly with x. 



Appendix E: Strong-coupling approach in the K x -model 

Here we give the details of the strong-coupling expansion 
around the limit where 1/4 of the plaquettes are decoupled 
from each other. Following the notation of Fig. 18, an isolated 
strong plaquette has the two triplet GS's with energy -3/16: 



) ay ® \t m )/3S , W S 



\t m ) ai ® \s) P s , (El) 



where \s) is the singlet and \t m ) are the three S=l triplets 
(to = ±1,0). So we have a GS manifold with degeneracy 
6 Np , where N p = N/4 is the number of the strong plaque- 
ttes. The splitting of this manifold by the inter-plaquette tC x 
terms can be captured by first order degenerate perturbation 
theory. The resulting effective model can be written in a con- 
venient form by introducing a pseudospin r = 1/2 and a spin 
T = 1 object in each strong plaquette. The direction of the 
pseudospin denotes which of the two types of triplets is taken, 
while the spin T = 1 object carries the physical spin of the 
plaquette. This approach follows closely in spirit the treat- 
ment by Lecheminant and Totsuka 51 of a very similar situ- 
ation where a different pair of plaquette triplets emerges at 
low energies in a two-leg ladder model. The inter-plaquette 
interactions are always of the type (S Q • S 7 ) (S^ ■ Ss) where 
a, /3, 7, S do not all belong to the same strong plaquette. In 
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FIG. 18: The plaquettized version of the square lattice. Thick diago- 
nal bonds denote the strong tC x -terms, while dashed ones are treated 
with first-order perturbation theory. 



particular there are two-plaquette interactions (see fig. 18) 

Vie = (s Q , ■ s 7 ) • s ff o = X) s ? s ^ s ™ s *' (E2) 

Vi y = (S a ■ Sy) (Sg ■ S$>) = y^^SpS™ S™Sg, ,(E3) 



and 4-plaquette interactions (see Fig. 18) 

V 2 = (Sa-SyHS^-S*™). 



(E4) 



So in the following we shall treat these interactions with first 
order perturbation theory, i.e. we shall write down an effective 
Hamiltonian which is formally given by 



H eS = PVP, 



(E5) 



where P — TJ P p is the projection operator that projects into 
the GS manifold, and 

p P = E (mm + \^T)(^T\) P ■ (£6) 

m=±l,0 

To continue we introduce a pseudospin r = 1/2, whose direc- 
tion specifies the type of the triplet, and a spin T = 1 object 
which carries the total physical spin of the plaquette. This 
defines the following mapping 

l^n-nt)®!!™), fofra*>-H+>®|T m >- (E7) 

It is straightforward to show that the pseudospin operators are 
the following scalar operators of the original spins of the pla- 
quette: 



T z = - ( — S a ■ S 7 + Sjj ■ Ss) 



(E8) 



t x = S a • Sp + S 7 • Ss = S Q • Ss + Sp ■ S 7 (E9) 
T y = Npg ■ S a x S 7 + N Q7 • S/3 x (E10) 
1 = -2(S a -S 7 + Sp-S s ) (Ell) 

where N Q7 = S a — S 7 and Nps = Sp — Ss- 

We now turn to the operators that are needed for the deriva- 
tion of the effective Hamiltonian. The following relations hold 
for single-site operators 

P S a P = P S 7 P =||>(|| (8) = [\ - r z ) 

PSpP = PS s P= |t)(t|®^T= (I + t z )®1t. 

For operators on the bond (/3, 7) we find: 

PiSZSy + SffS^P = \t x ® (Q nm - V m )(E12) 



P(Sp x S 7 )P = - Ty ® T 



(E13) 



where we have introduced the quadrupolar tensor of the spin 
T = 1 object: 



Qrirn _ rjvrirpm _|_ rpmrpn T{T -\- \^S nrn (E14) 

3 
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Similarly for the bond (a, 8): 



1. Symmetries 



P(S%S? + S™S£)P = X -t x <g> (Q nm - ^"" l )(E15) 



P(S a x S S )P = --T y ®T 



(E16) 



and for the bond (a, (3): 



P(S%Sp + S™S%)P = ~t x (8 (Q nm - ^ ,im (E17) 



P(S Q X S^P = ^Ty ® T 



(E18) 



Finally, the corresponding mappings for the bond (7, <5) are 
the same with that of (a, (3). 

Using the above mappings we may write down the first or- 
der effect from the inter-plaquette interactions. We have 

P V lx P = \t x t' x ® ( Q • Q' + I) +\T y T' y ®T-T' , (E19) 



where the vector Q contains the following elements of the 
quadrupolar tensor 30 

Q T = ( ±(Q xx -Q yy ),^Q zz ,Q xy ,Q yz ,Q xz ) ■ (E20) 
An equivalent way to write the above interaction is 

PV lx P = ^{t xT ' x - T y T y ) <g> n 5=0 + \{T X T' X + T yT ' y ) ® P l3 , 

(E21) 

where 

n s=0 = ^(Q-Q , + ^-T-T') ! (E22) 

3 

is the projector into the total singlet state, and 

Pij = \(Q ■ Q' + I + T • T') = -2 n s=1 + 1 , (E23) 

is the permutation operator which switches the states of two 
triplets. 30 

To switch from V\ x to V\ y , we need to map a <-> 7 in 
both plaquettes, which corresponds to 1 7) — > — | j) , and thus 
t± — ^ — r± and r x ,j/ — ^ ~ T x,y Since we need to do this 
in both plaquettes the minus signs cancel each other and thus 
P Vi y P has exactly thesame form as P V\ x P. 

The remaining 4-plaquette interaction term, denoted by V2 
above, gives 

PV 2 P = {\-T z ){\-T> z ){\+T>;){\+T>:') 



JL(T-T') (T" -T" 



(E24) 



It is interesting to note that the spin- 1 portion of this interac- 
tion looks exactly the same with the original four-spin inter- 
action tC x , but now the four sites have spin 1 and not 1/2. 



The original Hamiltonian has an SU(2)xSU(2) symmetry, 
i.e. we can rotate all spins labeled by a and 7 independently 
from the spins (i and i5. 65 The corresponding generators are 
the total spins of the two sublattices 



N p N p 

S J 4 = ^S Q7 (p), S B = ^Sps (p) 

p—1 p—1 



(E25) 



and the group elements are parametrized by two vectors fi^ 
and rtg: 

R(n A ,n B ) = R a (n A ) R b (n B ) = e -i(s A -n A+ s B -n B ) ; 

(E26) 

or, more exphcitly i?(r2 J 4, fl B ) — Yi P Ppi^A, ^b)> with 



R p ((l Ai n B ) = \t) P (t\®e~ lTp ' nB + \l) P (l\®e- 



(E27) 

As expected, all the effective terms written above retain the 
SU(2)xSU(2) invariance of the original /C^-model. To see 
this e.g. for the term P V\ x P, it is convenient to make use of 
Eq. (E21) rather than Eq. (E19). 



2. Quantum-mechanical solution for two neighboring 
plaquettes 

For two neighboring plaquettes pi = (1234) and p 2 = 
(5678) and for both open and periodic boundary conditions 
the QM GS is the following: 

IV') = i(| tita) - I UU)) ® \Ti = 1,T 2 = 1,T = 0) . 

(E28) 

The GS energy is -3/16, i.e. \ip) minimizes fully the interpla- 
quette interaction. We should note here that in this wavefunc- 
tion, the spin-1 objects make a singlet with (Ti • T2) = —2, 
and (Qi • Q2 + 2/3) = 4. The latter value is twice the one 
we would get by using a product state of two parallel director 
wavefunctions. 



3. Classical Variational solution to Vf = P(Vi x + V ly )P 

In the following we shall present a variational treatment of 
the terms V\ x and V\ y , disregarding P V2 P. The variational 
wavefunction is a product of plaquette wave-functions with a 
pseudospin portion and a spin-1 portion. The pseudospin-1/2 
portion is treated classically, i.e. it is parametrized by two 
direction angles 9, cf>. On the other hand, the spin-1 portion 
is parametrized by a complex vector d = u + ?v, with the 
constraints u 2 + v 2 = 1, and u • v = 0. 30 

A numerical minimization of this variational state delivers 
a GS in which the r = 1/2 pseudospins order AFM in a Neel 
state with their moments along the x-axis, while the spin-1 ob- 
jects develop a ferro-quadrupolar ordering without any dipolar 
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moment. Choosing the common director along the z-axis, we 
can write this state explicitly as 

N p 

|Z, Z> = ]J (l S )»^ l Z )^ + (-l) Q ' Rp |z)a7 ® I S )/m) I 
p=l 

(E29) 

where Q = (tt, 7r) is the Neel ordering wavevector of the 
pseudospins. We should note here that one can actually gen- 
erate a continuous family of equivalent variational wavefunc- 
tions by SU(2)xSU(2) rotations of |z, z). These states are of 
the general form 

\d a , d fa ) = Yl (|s) Q7 ® \d b )ffg + (-l) Q ' R »|d a ) Q7 ® \s)p S ) , 
p=i 

(E30) 

where now we have two (in general) different directors point- 
ing along d a and d(, in the two sublattices. By symmetry all 
these wavefunctions give the same variational energy and so 
we can choose any one of them to work with. 

In the following we take the state |z, z) and look at the ne- 
matic order parameters for the initial spin- 1/2 degrees of free- 
dom. Using the relations 

PS±S± P = \t)(t\®\(T ± ) 2 (E31) 

PS±S± p=\l)(l\®l(T ± ) 2 , (E32) 

and 

(t x =±- 2 \(- 2+ t z )\t x = ±-) = - (E33) 

(r^iilQ-^lr^i^i, (E34) 

and looking back at Fig. 18 we find: 

(S+S+) = (S+S+)?0, (E35) 

but e.g. (S^S+J = 0. There are two important points to 
note here. First, according to Eq. (E35), the above variational 
GS has a finite spin-nematic order parameter in both sublat- 
tices. This is however not in line with our ED results in the 
full JC X -model, since the low-E tower of states suggests that 
the spin-nematic ordering occurs only in one of the two sub- 
lattices (see Sec. XB). Secondly, the fact that (S'+S'+) is not 
equal to (S+ 3 ) is in line with the d-wave symmetry shown 
in Fig. 17. However, as we comment in Sec. X B, this result is 
to some extent artificial since the bonds (ai,7i) and (03, 74) 
become inequivalent once we plaquetize our lattice. 

Finally, we discuss the influence of an external magnetic 
field. Here the field reduces the symmetry of the Hamiltonian 
down to U(l)xU(l). Following Refs. [30,56,57], we expect 
that an infinitesimal field will pin the two directors on the xy- 
plane and will also induce a finite longitudinal dipolar moment 
in both sublattices. This picture is confirmed by a numerical 
minimization of the corresponding variational state in a field. 
In particular, we find that the spin-1 objects become fully po- 
larized at a critical field H c = K x /32. 



Appendix F: Tower of states spectroscopy 
1. Method 

Here we provide the group-theory derivation of the low- 
energy tower of states corresponding to some of the states that 
we encountered in this work. The method is quite general and 
is based on the following recipe from group theory. By apply- 
ing all elements of a symmetry group G on a given state |c) we 
generate a family (or orbit) of symmetry-equivalent states. By 
construction, this family of states provides a representation O 
of G and can be decomposed into IR's D" of G as follows 

0(5) =0 m a B a (g) (Fl) 

a 

where the number of times a given IR labeled by a appears in 
this decomposition is given by the well known formula 58-60 

1 1 sec 

Let us now assume that the starting state |c) is a classical mag- 
netic state (i.e. a site-factorized product of coherent states 
where each spin points to some fixed direction), and let us 
choose G = G r x SO(3) as the product of a real space group 
G r times the global SO(3) rotations in spin space. The gen- 
erated orbit O is then a continuous family of states and its 
decomposition into symmetry sectors can be obtained from 
Eq. (F2), if we can find a way to calculate the trace Tr[0(<?)] 
for each g € G. For finite systems this is a quite involved 
task given that one first needs to find an orthonormal basis 
that spans the orbit. One alternative method is to use projec- 
tion operators (which project onto specific rows of irreducible 
representations), but this relies on performing a number of nu- 
merical integrations over the group manifold for each symme- 
try sector and finite cluster separately. 

A much better alternative is to keep Eq. (F2), but make use 
of the thermodynamic limit. Here the generated orbit of the 
classical states themselves give an orthonormal basis because, 
in contrast to the finite-size case, two classical states that are 
related even by an infinitesimal global rotation are strictly or- 
thogonal to each other. One immediate consequence of this 
feature is that out of the continuous sum over g in Eq. (F2) 
only a much smaller (and usually discrete) subset of opera- 
tions give a non-vanishing trace. Specifically, these are the op- 
erations h that leave the state |c) invariant up to a global phase, 
namely h\c) — e l ^ h \c). These elements form a subgroup H c 
of G, called the stabilizer of |c). The above orthonormal basis 
is actually in one-to-one correspondence with the elements of 
the factor group G/H c , and Eq. (F2) reduces to 6667 

m « = T7T7 E ^ h X a {h) ■ (F3) 
1 c| heH c 

This is our central formula for the derivation of the tower of 
states for a magnetic state. 

A few important comments are in order here. First, the 
above derivation remains valid whether we take for G the 
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(a) collinear (b) orthogonal 

FIG. 19: Classical collinear and orthogonal states. 



TABLE III: Phase factors e i4>K that appear in the symmetry decom- 
position of the orthogonal and the collinear phase (see Fig. 19). Here 
h a stands for the spin-space portion of the given symmetry operation 
h of the stabilizer H c , 



orthogonal 


collinear 


h s 


e *<PK 


h s 




R(0) 

R x ,y(TT),Rz(±n/2) 


(-if 

(-1) N ' A 


R(0) 

R*{<t>) 

R*(it) 


1 
1 

(_l)iV/2 



complete group of the Hamiltonian or any subgroup of it. Sec- 
ond, one should always check whether the orbit generated by 
the elements of G covers the full family of classical states that 
we want to consider. For example, the family of the collinear 
states that are stabilized for large J2 / J\ comprises two sub- 
families of states. In the first (resp. second), the spins are 
FM (resp. AFM) along the line y' of Fig. 19 and AFM (resp. 
FM) along the line x' . These two sets are not related to each 
other by a global rotation in spin space nor by any primitive 
translation of the lattice. So if we take G =SO(3) alone or 
G = T® SO(3) (where T is the translation group) we gener- 
ate two distinct invariant orbits and each one must be decom- 
posed separately. By contrast, the above two sets of collinear 
states can be mapped to each other by a C4 spatial rotation. So 
if we include C4 in G, then the generated orbit will comprise 
both sets of states and so we have to do the decompositionjust 
once. 

Our final comment on Eq. (F3) is about the phase factors 
e *<Ph \y e ma y understand the origin and significance of these 
phase factors by taking as an example the 2-sublattice Neel 
AFM on the square lattice where half of the spins point along 
+z and the other half along z. Then one element of the 
stabilizer is a combination of a C-j spin rotation around e.g. 
the y-axis with a spatial operation that maps one sublattice to 
the other (e.g. a translation by one lattice constant). In this 
case we pick up a phase factor of e m l 2 from the action of 
the spin rotation on each spin site. With an even number 
of spins, this gives e l ^ h = (— 1) N / 2 . The same phase factor 
arises for a C2 spin rotation around any axis on the xy-plane. 
Overall we get a different decomposition for the square lattice 
Neel AFM depending on whether N = 4p or N = 4p + 2 
(where p is an integer). 68 

In the following we shall make use use of the central for- 
mula (F3) to decompose the orthogonal and the collinear fam- 
ilies of states into IR's of the group G = T C4® SO(3). Then 
we shall derive some further symmetry properties by includ- 
ing the glide operations (<7j|r) as well. 



2. Orthogonal phase 

For the orthogonal state we take the global spin orientations 
shown in Fig. 19(b). The stabilizer of this state consists of the 



following operations h 



{t 2 x n t 2 y m , tt n+1 tl m+1 } ®{E,C 2 } 



L x S 

f 2n+l ,2m ,2n .2m+ll 
^y ' ^y 

{tl n tl m , tf +1 t 2 y m+1 } ® {C 4 , C7 1 } ® R y (n) 
{tt n+1 t 2 y m , e t 2 y m+1 } ® {C4, C7 1 } ® R x (tt) 



where n,m = 1, ■ ■ ■ L/2, R xy z (n) are spin rotations by ir 
around the axes x, y, and z shown in Fig. 19. The phase fac- 
tors e l ^ h are provided in Table III. In the following we shall 
consider systems where the number of sites N is such that all 
phase factors reduce to 1 . 

Adding the contributions from the above operations h and 
using x S (R(0)) =25 + 1, and x 5 (tt) = (-l) s , we find that 
the only non-vanishing integers mk,;.s are those for k = 
and (tt, 7t), with 



m k =oA,s = \s lt{0 ,* } [25 + 1 + (-l) s (l + 2e* 



m-k=(Tr,7T),LS = T<^,{o,7r} 25 + 1 - (-1)" 



We note in particular that 



m c 



k,i 



"ikJ,s = 25+1 



(F4) 



so we have 2S+1 different total S-states for each given S. This 
important result can be also arise by a decomposition of the 
classical states into IR's of SO(3) alone. Since the orthogonal 
state breaks all global spin symmetries, and so the stabilizer 
consists only of the identity i?(0), whose character gives the 
(25 + 1) -multiplicity. 

We should remark here that the above symmetry struc- 
ture has also been reported for the orthogonal phase found 
inRef. [32]. 

Inclusion of glides (<7j|t).- So far we have not exploited 
the four non-symorphic operations (erj|r) of the lattice. The 
states with k = (tt, tt) belong necessarily to the 2-dimensional 
IR "(tt, 7r). {0, 7r}". So it is sufficient to discuss the states with 
k = only. In addition to the elements of the stabilizer listed 
above, now we may add the following symmetries that involve 
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TABLE IV: Group theory predictions for the content of the low-lying 
tower of states in the orthogonal phase. 



s 





1 


2 


3 


4 


0.A1 (I = 0) 


1 




1 




2 


0.A2 (I = 0) 






1 


1 


1 


0.B1 (I = tt) 




1 




1 


1 


0.B2 (I = tt) 






1 


1 


1 


(tt,tt).{1 = 0,tt} 




1 


1 


2 


2 


Total 


1 


3 


5 


7 


9 



a glide operation 

|£2n ^2m ^.2n+l f2m+l 



, C l ^*r^> ®(tr x ,v|r)(g)ii,(-7r/2) 

j.2n /.2m j,2n+l x2m+l 



ff\ tl m+1 } ® (<x XiS ,|t) ® ^(tt/2) 



where i? x / y '(7r) are spin rotations by tt around the x',y' di- 
rections shown in Fig. 19. By including these operations in 
the stabilizer group we get the following multiplicities at zero 
momentum: 



™k=o,r,s 



1 

To 



(2S+l)( X r (E)+x r (C2) 
+ (-l) S (x r (£)+X r (C 2 ) + 4 X r (C 4 )) 

+ 2(-l) s (x I > d |T)+ X I Wk)) 
+ 2 X S (7T/2)( X r (a x |T) + x r K|T) 



where T labels the IR's of C^ v shown in Table I. In Table V we 
show the specific multiplicities for the lowest four total spin 
sectors S. 



3. Collinear phase 

We turn now to the collinear phase with the global spin ori- 
entations shown in Fig. 19(a). As we mentioned above, by 
including the spatial Ca rotation in the group G we generate 
an orbit which includes both sets of collinear states. The sta- 
bilizer of any member of this orbit consists of the following 
operations 



{t 2 x n tf\ t 2 x n+1 t 2 y m+1 } ®{E,C 2 } ®{R x (m 



{t 



2n+l j2ro ^2n j2m+l 



} ®{E,C 2 } ®{R a (Tr)} 



where a is any axis perpendicular to the direction x of the 
spins (see Fig. 19), and (j> <G [0, 2tt). Again in the following we 
are considering system sizes N such that all phase factors e l ^ h 
corresponding to the collinear phase (see Table III) reduce to 



1. Wehave|{i? z (0)}| = \{R^)}\ = J t 



2tt 




2tt, and 



1 

27 



o 27T Jo 



27r sm(S 
dq> — 



= 1 
(F5) 



TABLE V: Group theory predictions for the content of the low-lying 
tower of states in the collinear phase. 



s 





1 


2 


3 


4 


0.A1 {I = 0) 
0.B2 (I = tt) 


1 
1 




1 
1 




1 
1 


(tt,tt).{1 = 0,tt} 




1 




1 




Total 


2 


2 


2 


2 


2 



for integer S (this is our case since we consider clusters with 
an even number of sites), and similarly 



1 

2Wo 



X S (i?a(7T)) 



1 

2W 



(-1) 5 = (-1)' 



(F6) 

So we find that the only non-vanishing integers mk.z,s ar e 
those with k = and (tt, tt) with 



1 



"lk=0,/,S = <5/,{0,7r}- 
Wk=(7r,7r),;,S = ^,{0,7r} 

and in particular 



1 -(-!)' 



i's = X] m k,;,s = 2 



(F7) 



k,/ 



so we have 2 different total S-states for each given S, one 
with I = and another with I = tt. This result also fol- 
lows from the decomposition of the classical states into IR's 
of SO(3) alone. Here the stabilizer consists of the U(l) ro- 
tations {R x (cf>)}, whose weighted integral over the character 
gives 1 . An extra factor of 2 comes from the fact that we have 
to consider each family of collinear states as a separate orbit 
(representation) since, as we mentioned above, they are not 
related to each other by SO(3) alone. 

Inclusion of glides ((X;|t).- Let us now exploit the four 
non-symorphic operations (<x;|t) of the lattice in order to 
make further predictions for the symmetry properties of the 
collinear tower of states under these operations. As we did 
above, it is sufficient to discuss the states with k = which 
appear for even values of S here. Since the collinear tower 
has only two states per spin sector S it is sufficient to look 
at the decomposition of the classical family of collinear states 
into IR's of Ci v alone, i.e. disregarding the operations in spin 
space and the translations. To this end one needs the corre- 
sponding (stabilizer) subgroup of C± v with the operations that 
leave a given collinear state invariant. It is easy to check that 
this is the C 2v subgroup which comprises in addition to E and 
C*2 the glides (od|r) and (ad' |t). Using the characters shown 
in Table I one then finds that only the IR's A\ and B 2 should 
appear in the tower of states when k = 0. This is because 
only A\ and B 2 remain invariant under the action of (c^t) 
and ((Trf'|x). So for each even value of S, we have one A\ 
and one B 2 state. The first is symmetric with respect to all 
non-symmorphic operations, while B 2 is even with respect to 
(f7<j|f ) and (ad 1 \t) but odd with respect to (a x \r) and (a y \r). 
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4. The spin-nematic variational state of App. E 3 

Here we derive the spin symmetry decomposition of the 
spin-nematic variational state |z, z) of Eq. (E29). To this end, 
we first need to find the subgroup of SU(2)xSU(2) rotations 
which leave the state |z, z) invariant up to a global phase. One 
set of operations that belong to the stabilizer is the U(l)xU(l) 
rotations around the z-axis 



So, the tower of states corresponding to the variational state 
|z, z) must consist only of sectors with even Sa + Sb- The 
numerical spectra for the /C^-model (see Fig. 16) show a dif- 
ferent symmetry pattern in the tower of states, so we believe 
that the variational state of Sec. E3 does not describe the GS 
of the /C^-model. 



R((t> a ,(j>b) = R a {4> a z)Rb((i>bz) 



(F8) 



Another set of operations is that of simultaneous it -rotations 
Coov x Coov in the two sublattices around any pair n a , rib of 
axes in the xy-plane 



i?(7rn a ,7rn;,) = i? a (7rn a )i? & (7rn b ) 



(F9) 



This can be shown by first noting that a 7r-rotation around 
an axis in the xy-plane reverses the sign of a triplet |i ), e.g. 
e -47rSl |to) = — |^o)- Hence, to restore a plaquette wavefunc- 
tion of the type 



1 



V2 



(|s)q 7 Ko)/35 ± \to) cfy\s) fit) 



(F10) 



(see Eq. (E29)), we must perform a ir -rotation in both sublat- 
tices A and B. Since this gives a minus sign for a single pla- 
quette, we get an overall phase of (— 1) p , where N p = N/4 
is the number of strong plaquettes. 

Having established the stabilizer of the state |z, z), we may 
now make use of Eq. (F3) in order to decompose it into sym- 
metry sectors of SU(2)xSU(2), which are labeled by the total 
spins Sa and Sb of the two sublatices. For clusters with even 
N/4, we get 



m s A ,s B = 2 ( 1+ 



Sa+Sb' 



5. U(l) spin-nematic phase 

The numerical tower of states shown of Fig. 16 consist of 
sectors (Sa, Sb) where one of the two spins is zero and the 
other is an even integer. The first suggests that the symme- 
try breaking occurs in one sublattice only, while the second 
suggests a spin-nematic state with U(l) symmetry. To demon- 
strate this we provide a simple derivation of the tower of states 
for a U(l) spin-nematic phase corresponding to a product of 
an even number of directors, all pointing along the z-axis. 

The stabilizer consists of global rotations around the z-axis 
as well as global 7r-rotations Coo V around any axis on the xy- 
plane. Given Eqs. (F5) and (F6) above for the corresponding 
characters, we arrive at 



1 



m s = , (F12) 

namely, only even sectors S participate in the tower. Further- 
more, the fact that we get a single state per even S in the tower 
of states is related to the U(l) symmetry of the spin-nematic 
state described here. If the directors were not collinear, as it 
happens e.g. in the antiferro-quadrupolar phase of the 5=1 
billinear-biquadratic model in the triangular lattice, 56 then a 
(Fl 1) more complicated multiplicity appears in the tower of states. 
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We should note here that by including the correction SEi (see 
above) from harmonic fluctuations around the classical up-down 
state we recover the correct energy — S(S + 1) J33 in the limit of 
x — ?> 0). Still however, the up-down state is very different from a 
singlet wavefunction. 

This follows by employing the theory of Ray representations or 
Herring's method (see e.g. Chapter 1 1 of Ref. [60] and references 
therein). 

This specific case is discussed in Sec. 1 1.2.5 of Ref. [59]. 
One may restore translational invariance by constructing a (zero- 
momentum) linear superposition of the four variational states cor- 
responding to the four different plaquettized versions of the square 
lattice. 

We should note here that the unperturbed part Ho has a much 
higher symmetry since we can do independent SU(2)xSU(2) ro- 
tations in every strong plaquette. 

In somewhat more technical terms, the classical state must belong 
to a given one-dimensional representation D x of H c (A is not al- 
ways the identity representation of H c because of the non-trivial 
phase factors e l ^ h ), and O is the so-called "induced representa- 
tion" D E j G. So what we are actually doing here amounts to a 
decomposition of D E f G into symmetry sectors of G. 
The derivation of Eq. (F3) also appears in App. B of Ref. [43] 
but without the phase factors e l< ^ h which were overlooked. This 
however does not affect the results of Ref. [43] since e l ^ h = 1 in 
all cases studied there. 

This signifies that there might exist different types of clusters that 
approach the thermodynamic limit in a different way. 



